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FOREWORD 


This report is one of several to be published from research con- 
ducted \ander NASA Contract NASS— 31773 entitled, "Relationships Between 
Severe Storms and Their Environment." This effort is sponsored by the 
NASA Office of Applications under the direction of Marshall Space Flight 
Center's Aerospace Environment Division. The results presented in this 
report represent only a portion of the total research effort. Data 
used in the report were taken from the AVE IV E3?periment conducted 
during the period beginning at 0000 GMT on 24 April 1975 and ending at 
1200 GMT on 25 April 1975. 
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ATMOSPHERIC ENERGETICS IN REGIONS OF INTENSE CONVECTIVE ACTIVITY 

By Henry E. Fuelberg^ 

Center for Applied Geosciences, Texas A&M University 
1. STATEMENT OF THE PROBLEM 

a. Introduction 

A knowledge of atmospheric energetics is basic to an understand- 
ing of the atmosphere. All atmospheric processes, ranging from the 
largest circulation systems down to the smallest turbulent motions, 
are accompanied by changes in energy. An understanding of the inter- 
action of the energy of these various scales of motion will help ex- 
plain and forecast many weather phenomena that presently are not un- 
derstood.' Thunderstorms are important daily weather producers which 
generate large quantities of energy, but the role of thunderstorms in 
the atmospheric energy budget is not understood well. This research 
assesses the impact of two large, intense squall lines on the synoptic- 
scale budgets of potential and kinetic energy during stages of squall 
line development. The use of rawinsonde data at 3- and 6-h intervals 
permits the temporal resolution of features that ordinarily would be 
unresolvable with the 12-h data that are routinely available. 

b. Previous studies 

The literature reveals that most research in atmospheric ener- 
getics has focused on large-scale systems, but suggests that small- 
scale systems, such as convective areas, may be important features 
of the large-scale energy budgets. 

1) The energetics of large-scale areas 

Oort (1964) , after reviewing the work of previous investiga- 
tors, described the atmosphere's primary energy cycle as beginning 
with the creation of mean available potential energy due to dia- 
batic processes, which is then converted to eddy available potential 
energy, eddy kinetic energy, and finally to mean kinetic energy. 

The annual cycle of atmospheric energetics in the northern hemi- 
sphere over a 5-year period has been described by Oort and 
Peixoto (1974) , and Peixoto and Oort (1974) . Numerous studies 
have considered average energy conditions over specific large 
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areas of the earth and/or for time periods shorter than five years; 
some of these are the available potential energy budget over the 
tropical Atlantic Ocean (Vincent, 1974), the kinetic energy balance 
over the tropical Pacific Ocean (Kung, 1974) , and the energetics 
over North America (Kung, 1970) . 

2) The energetics of limited areas 

Individual synoptic-scale systems with wavelengths of several 
thousand kilometers and lifetimes of several days are important 
components of the average energy budget computed over large areas 
and long time periods. Energetics research for such limited areas 
of the atmosphere is complicated by the presence of boundary-flux 
terms which vanish when the entire atmosphere is under consideration. 

Numerous studies dealing with the kinetic energy budgets of 
individual cyclones have been conducted. A summary of this work 
is provided by Kung and Smith (1974) who noted that a single 
mature cyclone system can account for 1/4 to 1/3 of the middle- 
latitude energetics. Recent studies of individual cyclones have 
been made by Ward and Smith (1976) and Kornegay and Vincent (1976) . 

A major conclusion reached by Vincent ^ (1974) in their 

kinetic energy study of Hurricane Celia was that the energy budget 
was dependent on the volume size used to contain the area of 
interest. This problem is overcome when values of terms in the 
energy budget equations are presented throughout the spatial and 
temporal extent of a cyclone instead of presenting them in terms 
of averages only. Such results have been given by West (1973) and 
Ward and Smith (1976) v/ho could track in both time and space the 
important energetic processes of developing cyclones. Smith (1973a) 
studied the kinetic energy budget of an anticyclone and found it 
to be a tranquil circulation system when compared to cyclones. 

Several studies have considered various aspects of the poten- 
tial energy budget of limited areas. Smith and Horn (1969) 
described the complete available potential energy budget for North 
America during March 1962, while West (1973) studied the budget of 
total potential energy during two cases of cyclogenesis. Goinmel 
(1973) presented values of terms describing the conversion of 
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available potential energy in Hurricane Celia, while Kung (1974) 
has done the same for an area over the tropical Pacific Ocean. 

The effect of geographical heat sources and sinks on the generation 
of available potential energy has been studied by such authors 
as Gall and Johnson (1971), Bullock and Johnson (1971), and Anthes 
and Johnson (1968) . 

3) The energetics of areas of convection 

Research indicates that areas of intense convection produce 
effects on the surrounding atmospheric volumes that can be detected 
using synoptic-scale data. Axobert (1957) showed that latent 
heat release associated with thunderstorms produced increases in 
large-scale horizontal convergence below the level of maximum 
condensation and enhanced horizontal divergence above that level. 

He found decreases in values of geopotential height in the lower 
troposphere and increases in the upper troposphere. Ninomiya 
(1971a and b) , using satellite pictures, found that the pre- 
existing flow at the cirrus level over tornado-producing thunder- 
storms changed into outflow as the thunderstorms developed. The 
existence of a midtropospheric warm core and a significant field of 
convergence below 700 mb was observed using synoptic- scale rawinsonde 
data, and the downward convective transport of horizontal momentum 
was shown to intensify the low-level jet stream. Each of these 
phenomena was attributed to latent heat release. Fankhauser (1971) 
indicated that mature thunderstorms may divert and distort mid- 
tropospheric air motion in a manner similar to that of solid ob- 
jects. While several authors have indicated that latent heat 
release is important in offsetting the frictional effects in a 
developing cyclone (Bullock and Johnson, 1971; Danard, 1964, 1966), 
Tracton (1973) and VJest (1973) have suggested that latent heat 
release plays a critical role in initiating cyclogenesis in some 
instances - 

Since research has shown that areas of convection affect 
synoptic-scale fields of kinematic and thermodynamic quantities, 
the energy budget should be affected also, but surprisingly little 
work has been conducted to describe the energy changes. Danard 
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(1964, 1966) used mamerical methods with and without the inclusion 
of latent heat to compute the effects of a large precipitation area 
on certain energy budget terms. The production of kinetic energy 
due to rising air and cross-contour flow was enhanced considerably 
by latent heat release whose effect was of the same order of magni- 
tude as that caused solely by dry adiabatic processes. The kinetic 
energy budget of atmospheric volumes surrounding trade wind cloud 
clusters in the Pacific Ocean has been investigated by Williams 
(1970) who found them to be vertical transporters of kinetic energy 
with little internal generation or dissipation. Kinetic energy was 
found to be imported in the lower half of the troposphere, trans- 
ported upward within the convective clouds, and exported in the 
upper tropospheric outflow. The kinetic energy budget for clear 
areas was considerably different from that of the clusters. 

The kinetic energy budget of volumes enclosing areas of con- 
vection within the central United States has been presented by 
Mclnnis and Kung (1972) using data from the rawinsonde network 
operated by the National Severe Storms Laboratory (NSSL) . Their 
results were reported in terms of spatial averages for the entire 
NSSL area at various pressure levels and times. Kinetic energy 
generation and dissipation were found to be maxima in the lower 
troposphere and at the jet stream level. Horizontal inflow and 
vertical outflow occurred in the lower and middle troposphere while 
horizontal outflow and vertical inflow were found near the jet 
stream level. 

Much remains to be learned about the energetics of all scales 
of motion, but additional studies dealing with scales smaller 
than synoptic, such as convection, would appear to be especially 
useful. 

c. Objectives 

The objective of this research was to test the hypothesis 
that systematic changes occur in the synoptic-scale energy budgets 
of atmospheric volumes enclosing areas of intense convection, and 
that these changes are related to the life cycle of the convective 
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areas and can be tracked in space and time. 

To accomplish the objective it was necessary to: 

1) Use energy budget equations that define energetics in a 
limited region of the atmosphere and that allow the computation 
of boiindary fluxes of energy and the interaction of the convective 
area with its surroundings. 

2) Develop analytical procedures for evaluating the energy 
equations that yield accurate results without requiring excessive 
smoothing that would destroy significant features. 

3) Compute spatial fields of terms in the energy budget equa- 
tions at various pressure levels during the formation, maintenance, 
and dissipation of areas of intense convection. This was done to 
relate the locations and magnitudes of the energetic processes to 
the locations and life cycles of the convective activity and other 
weather systems. 

4) Compute average energy budgets for volumes of various 
sizes that enclose areas of convection and nonconvection so that 
results of this investigation could be compared with those of 
meteorological systems of large scales. The influence of intense 
convection on these systems could then be assessed. 
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2. THEORETICAL DEVELOPMENT 


The energy budget equations used in this study were proposed by 
West (1973) based on suggestions by Smith (1970). Individual terms 
of the equations represent physical processes that lead to changes 
in atmospheric energy and provide a means to study the interaction 
of the limited volumes \inder consideration with larger voltme.s. 
Highlights of West's (1973) derivation of these equations along 
with certain modifications are given in this section. 

A limited atmospheric volume, V, possesses internal energy, I, 
due to its temperature, T; gravitational potential energy, P, due 
to its height, z; and kinetic energy, K, due to its vector velocity, 
V, These forms of energy are given by 



( 1 ) 

C2) 

C3) 


where c^ is specific heat at constant volume, p is density, g is the 
acceleration of gravity, and V" is a dummy variable of integration 
for volume integrals. 

Wiin-Nielsen (1968) , using the equation of mass continuity, 
presented the following equation that is applicable for a volume 
extending over a limited region of the atmosphere: 
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where b is any arbitrary scalar defined in the atmosphere, t is time, 
and ^ is the three-dimensional del operator. This equation may also 
be written as 


f f ~f ^*(5=PV)dV" + f 

-'v Jv Jv Jv 


pvVbdV", 


(5) 


or 


-f 

V V Jv 


V* (pbv)dv". 


( 6 ) 


V V 

The. last term in (6) is a boundary term that would vanish if inte- 
gration were performed over the entire atmosphere. 

When (6) is applied to (1) and the boundary of the volume is 
fixed in space, the time rate-of-change of internal energy is given 


by 


9l 

9t 


-'V Jv 


(pc^T^)dV". 


(7) 


dT 


An expression for c^ can be obtained by combining the thermo- 


dynamic equation with the continuity equation, 
then be written as 


Equation (7) can 


9t 

dH 


= f ^ - r 

J V 


VdV" - 


V-(pc TV)dV", 

V 


( 8 ) 


where ^ is diabatic heating per unit mass and p is atmospheric 
at 

pressure. 

When (6) is applied to (2), the time rate-of-change of 
potential energy for a fixed volume is given by 


9t 


I gpwdV" - I 


(gpzV)dV**, 


(9) 


where w is the vertical component of V. 

Finally, when (6) is applied to (3) and substitutions are made 
from the vector equation of motion, the time rate-of-change of 
kinetic energy for a fixed volume can be written as 
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3k 

9t 



V*^pdV” 



pgwdV" 



( 10 ) 


where F is the three-dimensional frictional force. 

The first term on the right hand side of (10) and the second 
term on the right hand side of (8) represent the conversion be- 
tween I and K in a limited volume if there is no boundary effect. 
These two "conversion" terms and the boundary term are related by 



r 

+ I V*VpdV" = 



( 11 ) 


The term on the right hand side of (11) usually is referred to as 
the pressure work term or the boundary work term and vanishes 
when integration is performed over the entire atmosphere if the 
lower boundary is rigid. Smith (1970) has shown that the boundary 
term also can be interpreted as the interaction of a limited volume 
with the surrounding volume. When (11) is substituted into (8) , 
one obtains 


8l 

8t 




V-VpdV" 


j ^*(pc TV)dV". 
V 


J ^*pVdV" 
V 


( 12 ) 


The hydrostatic approximation is applied since it already has 
been used in computing the rawinsonde data that serve as input 
for the budget equations. The assumption is a potential source of 
distortion in this research since convective cells undergo non- 
hydrostatic vertical accelerations, but steps were taken to lessen 
this potential problem (see Section 3) . 

The vertical component of K becomes irrelevant in the energy 
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budgets when the hydrostatic approximation is made and may be ne- 
glected (Wi in— Nielsen, 1968) such that the symbol K now represents 
only the horizontal component of kinetic energy. The vertical ve- 
locity required to maintain hydrostatic equilibrium is then denoted 


by Wj, 

9l 

9t 


With this assumption the equations for I, P, and K become 

8p 


9t 


9k 

9t 


- 

( 


IfdV + 

1 + 


J V 

V*pVdV" - 

1 V* (pc TV)dV", 


fw ^ 

J/r 9 . 


dV" 


V 


pgw dV" 
R 


I 


■ - I,': 


_pdV" 


^•(gpzV)dV", ^nd 


1 V •F^dV" 

J V 


(13) 


(14) 



where 


v^*v 
2 2 ^ 
P — ;r- V 


K = 


->■ 

I 

-^v 


dV", (15) 


(16) 


dV", 


is the two-dimensional del operator, F is the two-dimensional 
frictional force, and is the horizontal wind vector. The poten- 
tial and internal energies for the entire vertical column of a 


hydrostatic atmosphere remain in a constant ratio (Petterssen, 

1956) and are given by P = — I where R is the gas constant. 

Cv 

Total potential energy, JI, is defined as the sum of internal 
energy and gravitational potential energy and is given by 


n = 


X' 


PTC dV". 
P 


(17) 


When (13) and (14) are added, one obtains the time rate-of-change 
of total potential energy for a fixed volume. 


an 

3t 


X 

- Ij 


p g dV- 


• (pc TV)dV" - 

V 


^ I 


X 


^•pVdV" 


V* (pgzV)dV". 


(18) 



The sum of the last two terms on the right hand side of (18) was 
approximated by the term JT (pc TV)dV", where c is specific heat 

p p 

at constant pressure. Sample calculations indicate that values 
obtained using the last two terms of (18) are smaller than those 
obtained using the approximation in the lower and middle troposphere 
and are larger than the approximate values in the upper troposphere 
and stratosphere. However, values obtained from the two methods 
generally differ less than 25% so that conclusions based on the 
approximation are not affected significantly. 

Transformation into the x, y, p, t system provides the final 
equations for the horizontal component of kinetic energy and total 


potential energy of a fixed volume which are, respectively. 


and 


■■ill 


2 2 


dpdA", 


n 


c 

g 


II 


T’dpdA", 


(19) 


(20) 


'A-p 

where A is the horizontal cross-sectional area of V and A'* is a 

dummy variable of integration corresponding to A. 

When (18) with its approximation and (15) are transformed into 

9 ( Db V o ) 

the X, y, p, t system, terms of the type bV^*V p, — ^ ‘V z, and 

3 2 p 9z p 

boo occur, where b is an arbitrary scalar. These terms can be 

neglected because they are generally two to four orders of magni- 
tude smaller than the remaining terms. Equations (18) and (15) 


then become, respectively, 
G 


‘'K 


-■Hi 


dH 

dt 


dpdA” 


PW 


Jl 


zdpdA" 
2 p ^ 



3p' 


) 


(T03) I dpdA" 


VF 


il 


c 

g 



^ • 


TV dpdA" 


C 

_ _£ 

g 
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(TO)) dpdA" 


( 21 ) 



11 


and 


3t 



K 


ladpdA" 

2 p 


D 


r 

+ 



1 -y 

— V^-F^dpdA" 
p 2 2 


^ 


( 22 ) 



where (u = ~ and is the two-dimensional del operator on a constant 
dt p 

pressure surface. 

The terms in (21) and (22) represent the following physical 

processes. The terms are denoted by symbolic notation which does 

not include negative signs. Term G in (21) represents generation 

of potential energy due to diabatic effects. Conversion between 

kinetic and total potential energy due to crosg-contour flow is 

denoted by term C in (21) and (22) . Term PW of (21) defines 
K 

pressure work at the boundary which also represents interaction 

of the limited volume with the surrounding volume. The horizontal 

and vertical components of the boundary flux of total potential 

9IT 

energy are given by terms HF^ and VF^ of (21) . Finally, is the 
local change of total potential energy that results from the pro- 
cesses just described. Sairple calculations indicated that term 
9II 

generally was several orders of magnitude smaller than the 
o t 

largest terms in (21) . It will therefore not be included in siJb- 
sequent discussions. 

Local changes in the horizontal component of kinetic energy, 

9k 

in (22) are due to four processes. Term C has been described 
d t K 

previously as conversion between the two forms of energy. The dis- 
sipation of kinetic energy is denoted by term D. The exact in- 
terpretation of this term is dependent on its method of computation 
and is described in greater detail in Section 3. Terms HF and VF 
represent the horizontal and vertical boundary fluxes of kinetic 


energy. 
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Another quantity that is of interest is the release of available 
potential energy which is given by 



where a is specific volume. This term, which Smith (1969) calls 

A-conversion, represents changes in available potential energy 

due to the vertical redistribution of mass within the volume. 

Many pxiblished papers have used this term to represent energy 

conversion instead of the te3rm C , which Smith (1969) calls 

K. 

K- conversion. Term represents direct conversion between 
potential and kinetic energy only if there is no boundary effect. 
The interaction between the limited volume and the surrounding 
volume can be inferred by comparing the values of terms and 

Figure 1 is a schematic representation of the energy processes 
that are described in (21) and (22) ; it is a modified version of a 
figure presented by Smith (1970) . The solid boxes represent the 
limited atmospheric volume while the dashed boxes denote the 
surrounding volume . 

Equations (19) to (22) are the energy budget equations that 
were used to investigate the relationship between energy processes 
and areas of intense convection. 
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3. DATA AND ANALYTICAL PROCEDURES 


AVE IV experiment 

A data set at time intervals smaller than the usual 12-h inter- 
val and encompassing a period of intense convective activity was 
needed for this study. Data from the fourth Atmospheric Variability 
Experiment (AVE IV) conducted on 24-25 April 1975, and sponsored by 
the National Aeronautics and Space Administration (NASA) met these 
requirements . 

1) Rawinsonde data 

Forty- two rawinsonde stations participated in the AVE IV ex- 
periment; these are shown in Fig. 2 and listed in the Appendix. 
Soundings were taken at nine times — 24 April at 0000 GMT, 0600 GMT, 
1200 GMT, 1500 GMT, 1800 GMT, and 2100 .GMT, and on 25 April at 0000 
GMT, 0600 GMT, and 1200 GMT. 



Fig. 2. 


Rawinsonde stations participating in the AVE IV 
experiment. 
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Data reduction procedures used to process the AVE IV rawin- 
sonde data were designed to obtain the most accurate results 
possible (Fuelberg, 1974) , The raw angle and ordinate data were 
checked for errors prior to calculating the soundings, and computed 
soundings then were rechecked with corrections made as required. 

Data were given at 25-mb intervals by Fucik and Turner (1975) . 
Estimates of RMS errors of the thermodynamic quantities are given 
in Table 1 (Scoggins and Smith, 1973; Fuelberg, 1974). 

Table 1. RMS errors of thermodynamic (quantities for the AVE IV data. 


Parameter 


Approximate RMS error 


Temperature 1®C 

Pressure 1.3 mb surface to 400 mb 

1.1 mb between 400 and 100 mb 
0.7 mb between 100 and 10 mb 

Humidity 10% 

Pressure Altitude 10 gpm at 500 mb 

20 gpm at 300 mb 
50 gpm at 50 mb 


Recent studies have indicated that errors in the thermodynamic quan- 
tities may be even smaller than those given in Table 1 (Lenhard, 

1973; Brousaides, 1975) . 

An error analysis conducted by Fuelberg (1974) gave RMS errors 
of scalar wind speed and wind direction for the AVE IV baroswitch 
contact data. Table 2 presents these RMS errors for elevation angles 
of 40° and 20°. RMS errors for the smoothed 25-mb data that were 
used in this study would be somewhat smaller. 

Original strip charts from all rawinsonde soundings were checked 
carefully to determine - if sondes entered thunderstorm updrafts or 
downdrafts. Data for four soundings 'were removed from the original 
data set because the sondes apparently had entered violent up- 
drafts. The potential for distorted results due to nonhydrostatic 
accelerations was thereby considerably reduced. 
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Table- 2. RMS errors of baroswitch contact wind data for the AVE IV 
data at elevation angles of 40° and 20°. 


Elevation Angle Elevation Angle 

Level 40° . 20° 40° 20° 

RMS Direction Error RMS Speed Error 


700 

mb 

1.8° 

U) 

00 

o 

0.5 

m 

-1 

s 

1.0 

m 

-1 

s 

500 

mb 

2.5° 

5.6° 

0.8 

m 

-1 

s 

2.0 

m 

-1 

s 

300 

mb 

3.1° 

7.5° 

1.0 

m 

-1 

s 

3.8 

m 

-1 

s 

100 

mb 

6.2° 

15.0° 

2.0 

m 

-1 

s 

5.7 

m 

-1 

s 


2) Surface data 

All available surface data for the AVE IV experiment were ob- 
tained from the National Climatic Center. Figure 3 shows locations 
of the 310 surface stations used. Vector wind, temperature, dew 
point temperature, and surface pressure were read, keypunched, and 
checked carefully for errors. 



Fig. 3. 


Locations of surface stations for the AVE IV area. 
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3) Digital radar data 

Manually Digitized Radar (MDR) data were obtained from NOAA's 
Techniques Development Laboratory to determine accurately the in- 
tensity and position of the radar-observed convection during the 
AVE IV experiment. The MDR grid network is shown in Fig. 4. MDR 
data are coded with a single digit from 0 to 9 to indicate areal 
coverage and echo intensity within blocks that are approximately 
83 km on a side. Table 3 is an explanation of the MDR code given 
by Foster and Reap (1973) . 


Table 3. Manually digitized radar (MDR) data code (Foster and 
Reap, 1973) . 


Code No. 

Maximum 
Observed 
Vipl Values 


Coverage 
In Box 

Maximum 
P^ainfall 
Rate (in h“l) 

Intensity 

Category 

0 

No Echoes 





1 

1 


Any VIPl 

<0.1 

Weak 

2 

2 

< 

50% of VIP2 

0.1-0. 5 

Moderate 

3 

2 

> 

50% of VIP2 

0. 5-1.0 

Moderate 

4 

3 

< 

50% of VIP3 

1. 0-2.0 

Strong 

5 

3 

> 

50% of VIP3 

1. 0-2.0 

Strong 

6 

4 


50% of VIP3 
and 4 

1. 0-2.0 

Very Strong 

7 

4 

> 

50% of VIP3 
and 4 

1. 0-2.0 

Very Strong 

8 

5 or 6 


50% of VIP3, 
4 , 5 , and 6 

>2.0 

Intense or 
Extreme 

9 

5 or 6 

> 

50% of VIP3, 
4, 5, and 6 

>2.0 

Intense or 
Extreme 


^Video Integrator Processor 


Plots of the MDR data were made each hour for 3-h periods cen- 
tered on each of the nine rawinsonde observation times. The three 
plots were then combined into a single chart for each of the nine 
times by using the highest coded value reported for each block. 


L 



Fig. 4 


Manually Digitized Radar (MDR) grid network 
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b. Analytical procedures 

Analytical procedures can determine the success or failure of . 
an experiment and, therefore, must be considered carefully. This 
section describes the procedures that were judged most advantageous 
for this study and briefly discusses error estimates of terms in the 
budget equations using these procedures. 

1) Objective analysis and smoothing 

The computation of terms in the energy budget equations is 
simplified if data are interpolated from randomly- spaced rawin sonde 
stations to equally- spaced grid points. The grid network used (Fig. 5) 
is centered over the AVE IV data area and has a spacing of 158 km. 

Barr et al . (1971) have shown theoretically that a 169-km grid in- 
terval incorporates as much detail as is justified from the rawin- 
sonde network over the United States . 

Hirnian analysis of data is still preferable to any objective 
scheme that has been devised, but the number of analyses required 
rendered hand analysis unfeasible for this study. There were sev- 
eral considerations in the choice of an objective analysis and 
smoothing procedure from the several that are described in the 
literature. Such a procedure should interpolate data accurately 
from stations to the grid without creating fictitious waves or 
destroying real waves present in the data. Also, the desired 
results of the procedure should be considered carefully with respect 
to the input data density. Since the average spacing of rawin- 
sonde stations over the United States is about 400 km, features 
with wavelengths shorter than 800 km can never be described com- 
pletely by any objective analysis scheme. Barr ejt (1971) con- 

cluded that consistent qualitative evaluation of features with 
wavelengths less than 1700 km is unlikely. 

An objective analysis scheme by Barnes (1964) was used in this 
investigation. The procedure is commonly referred to as successive 
corrections to a first-guess field. Data from each rawinsonde 
station were allowed to influence grid points within a scan radius 
of three grid distances while four iterations were allowed. 

To suppress small waves which cannot be tracked consistently. 
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Fig. 5. Grid used for. numerical computations. 


as well as those which might arise due to the analysis scheme and 
to reduce random errors, a nine-point filter by Shuman (1957) was 
applied to the analyzed fields. The final result retained approxi- 
mately 90% of the anplitudes of wavelengths of 1600 km, while re- 
tention of wavelengths equal to 800 km and 400 km was approximately 
63% and 17%, respectively. These fields appear to retain as much 
detail as can be justified from the input data and have shown very 
good agreement with hand analyses . 

The surface data were analyzed using a scan radius of 2 grid 
units with four iterations. This produced smooth fields from the 
high density surface data that meshed with the larger- scale rawin- 
sonde data. 
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Gridded analyses of height, temperature, wind components, and 
moisture content were produced at 18 levels, i.e. the surface and 
at 50-mb intervals from 900 mb to 100 nto, for each of the nine time 
periods. Winds at the 18 levels were averages over 50 mb to reduce 
further random errors. These gridded values were stored on a com- 
puter disk and formed the working data set for all phases of the 
research. 

2) Numerical evaluation of equations 

i. . Finite differencing and integration . Centered finite dif- 
ferences were used to compute horizontal derivatives, and all ver- 
tical derivatives except those at the surface where forward differ- 
ences were used and at 100 nib where backward differences were used. 
Time derivatives also were evaluated using centered differences 
where possible, but forward differences at the first period and 
backward differences at the last period were required. 

Terms in the energy budget equations (19) to (22) were com- 
puted at each of the 18 levels previously described. These re- 
sults then were integrated over 50-mb layers from the surface to 100 
mb using the trapezoidal rule. Integrated results were stored on a 
computer disk for future use. Printed spatial fields of the energy 
budget terms were made for the layers surface to 700 mb, 700 to 
400 mb, 400 to 100 mb, and surface to 100 mb by integrating over 
these layers. Numerical evaluations were performed on the Amdahl 
470/V6 computer at Texas A&M University. 

Error estimates . Random errors in input data lead to 
errors in computed terms of the energy budget equations. Estimates 
of RMS errors expected in the results were made using the propagation 
of error method (Young, 1962) . 

The propagation of error method is briefly derived as follows. 

If F is a function of X, Y, and Z, and if these are in error by 
Ax, Ay, and Az, then F is in error by an amount Af which is given 
by the Taylor's series as: 



+ 


dF 

9Y 


Ay 



Af 


(23) 
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The equation is valid if Ax, Ay, and Az are small compared with 
3 f 8f 

and it is asstamed that all terms involving squares, higher 
powers, and cross products of Ax, Ay, and Az are negligible, and 
that higher derivatives are small in value. After squaring both 
sides of (23) and replacing each term by its average value, the 
following equation is obtained when one assumes that errors in X, 

Y, and Z are independent and normally distributed: 

where a denotes RMS error. 

Although the best estimates of RMS errors for the rawinsonde 
contact data were used to compute error estimates of the energy 
budget terms, several considerations suggest that the results prob- 
ably should be used only in terms of order of magnitude. First, 

RMS errors for objectively analyzed and smoothed fields of data are 
difficult to determine. Secondly, the effect of finite difference 
and integration procedures on errors is difficult to describe. 
Finally, the extent to which the input data agree with the assump- 
tions of the propagation of error technique is not known fully. 
Valuable insight into the error problem was obtained in spite of 
the difficulties involved in error estimation. 

The RMS error of each term in the kinetic energy equation, (22) , 

was compared with the range of values expected for the term. The 

RMS errors of all terms were found to be less than the e3q>ected 

values in the lower and middle troposphere. The cross-contour flow 

term, C , is the most sensitive at all levels, and its error ap- 
K 

proaches the expected values near the jet stream level. The er- 
rors of the remaining terms of the kinetic energy equation near 
the jet stream level are small compared to expected values. 

Recent results by Vincent and Chang (1975) and Kornegay and 
Vincent (1976) indicated that random errors in input data do not 
affect seriously the results of kinetic energy studies. After 
kinetic energy budgets were computed using real data, random errors 
approaching the extremes to be expected were introduced into the 
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wind and height data and the budgets were recomputed. Sign uncer- 
tainty was found to be greatest when the original result was near 
zero, but percentage changes between real and erroneous results de- 
creased as the original values become larger. Values of individual 
terms using erroneous data rarely exceeded original values by a fac- 
tor of two. Finally, the trends of results and their interpreta- 
tions almost always remained the same regardless of the errors intro- 
duced. These findings suggest that results based on the carefully 
checked AVE data will be valid. 

Random error estimates of most terms in the potential energy 
budget equation, (21) , are relatively small near the surface when 
compared to expected values, but increase with altitude. In the 
upper atmosphere, terms HF^^, VFj^, C^, and possess relatively 
large errors compared to the expected values of the terms, but based 
on the results of Vincent and Chang (1975) and Kornegay and Vincent 
(1976), no problems in interpreting the results should occur. The 
RMS error of the pressure work term, PW, in (21) is about two orders 
of magnitude larger than its expected value. The problem occurs 
because this term is the s\m of two component terms which are of 
nearly equal magnitude but opposite sign. Great caution must be 
used in reaching conclusions based on term PW. 

The difficulty in estimating errors for this type of energy study 
again should be emphasized. Mclnnis and Kung (1972) stated that 
judgements on validity of results usually are based on the systematic 
nature of those results and comparisons with other investigators. 

Kung and Smith (1974) pointed out that the effect of random errors 
can be reduced by averaging results over a large number of points. 
Further comments on error and interpretation of results will be made 
in later sections. 

iii. Computation of vertical motion . Large-scale vertical 
motion cannot be measured directly, and no method currently available 
for determining vertical motion is completely accurate. The kine- 
matic method was used in this research because it involved the least 
stringent assumptions and produces good results. Details of the 
procedures were given by Wilson and Scoggins (1976) . Terrain- induced 
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vertical motion was included, and a correction scheme by O'Brien (1970) 
was applied so the values of vertical motion at 100 mb would equal the 
values obtained by the adiabatic method. Wilson and Scoggins adjusted 
values of vertical motion to zero at 100 mb. The adjustment factor 
significantly affects vertical motion in levels above about 500 mb. 
Application of an adjustment factor is necessary because the accuracy 
of wind data and resulting divergence calculations decrease with al- 
titude. Adiabatic values at 100 mb were chosen because they are ob- 
tained independently of kinematic values and are more realistic than 
an assumption of zero at each grid point. 

The kinematic method has been used widely in previous research. 
Vincent ^ al . (1976) suggested that it is better than the quasi- 
geostrophic form of the omega equation. Further support for the 
kinematic method has been given by such investigators as Chien and 
Smith (1973), Smith (1971), Fankhauser (1969), and Kung (1973). 

An unpublished study by Wilson at Texas A&M University indicated 
that values of kinematic vertical velocity relate better to areas 
of precipitation during the AVE IV e:q)eriment than do values of 
adiabatic vertical velocity. 

iv. Frictional dissipation of kinetic energy . Term D in (22) 
was determined as a residual. Since dissipation of kinetic energy 
in large-scale circulations is not understood fully (Smith and 
Adhikary, 1974) , computation as a residual is the only objective 
way to compute it without employing specific theories regarding 
its nature (Kung and Smith, 1974) . While D conceptually represents 
thermodynamic frictional processes, Kung and Smith pointed out that 
when computed as a residual, it represents the sink of energy from 
the larger grid scales of motion which should serve as input to the 
smaller subgrid scales. It also contains errors accumulated from 
other terms. Computation of dissipation as a residual has been used 
by Mclnnis and Kung (1972) in their kinetic energy study of convec- 
tive areas and by most other investigators of energetics, 

V, Diabatic heating . Diabatic heating, term G in (21) , is 
composed of latent heat release, sensible heat transfer, and ra- 
diational effects. The literature suggests procedures for estimating 
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each of these components ^ but each procedure requires specific as- 
sumptions which may or may not be true at a particular point in space 
and time. Term G was evaluated from the thermodynamic equation 
<3H _ 9T . _ ^ 


dt 


= c 


p dt 


p 2 


(si - 


(25) 


which gives the net effect of the various components of diabatic 
heating. This procedure has been used by Smith and Horn (1969) and 


West (1973) . 



26 


4 . SYNOPTIC CONDITIONS 


Although no outstanding pressure or frontal activity occurred 
during the AVE IV experiment, two large areas of intense convective 
activity made the period ideally suited for this study. The synoptic 
analyses used in the following discussion are taken from Fucik 
and Turner (1975) and are similar to facsimile analyses made by the 
National Weather Service. 

Two air mass types were distinguishable at the surface over the 
AVE IV network at the beginning of the experiment at 0000 GMT on 24 
April (Fig. 6) . Warm, moist, maritime tropical (mT) air was moving 
northward over the network with strong southerly flow around a high 
pressure center located about 500 km off the coast of the Carolinas. 
This air mass covered almost two-thirds of the network extending 
from central Texas and Oklahoma eastward through all of the Gulf 
Coast and Middle Atlantic States and northeastward into the Ohio 
Valley. 

A cold front, extending southwestward into Kansas from a 
moderately strong cyclone (1000 mb center) over northern Michigan, 
separated the mT air from the cooler and drier continental polar 
(cP) air moving southward over the Northern Plains States. A warm 
front, extending southeastward from the cyclone into Pennsylvania, 
separated cP air over New England from mT air in the Ohio Valley. 

A second cyclone (1000 mb center) was located in central Kansas with 
a cold front extending southward into west Texas. Surface temper- 
atures over the AVE IV network ranged from about 25°C along the 
Gulf Coast to about 10 °C in the Northern Plains States, and dew 
point temperatures ranged from about 20°C to about -10°C. 

The middle and upper tropospheric flow pattern was basically 
zonal during the entire experiment with the polar jet stream ex- 
tending roughly west-east from the Northern Plains States to New 
England and the subtropical jet stream being located along the 
Gulf Coast States. TWo clearly distinguishable short wave per- 
turbations with wavelengths of about 1500 km moved through the 
zonal flow during the experiment period and were responsible for 



16 


Fig. 6 
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Synoptic charts for 0000 GMT, 24 April 1975 
(Fucik and Turner, 1975) . 
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much of the significant weather. There were also numerous shorter 
waves (750-1000 km) which could be resolved only partially by the 
data network. The first well-defined short wave was associated 
with the cyclone over northern Michigan at the beginning of the 
period. The second wave was still west of the AVE IV network at 
that time. 

The first major outbreak of severe weather began just prior to 
the start of the AVE IV experiment and was associated with the cy- 
clone in central Kansas and the warm front extending northeastward 
into northern Missouri. Maps of MDR data are given in Fig. 7 for 
the nine time periods. The area grew rapidly in strength and in- 
tensity so that tornado watches were issued for thh area at 0000 
GMT on 24 April. Radar tops extended as high as 18.3 km (60,000 ft) 
with MDR values as high as 9 in this unstable area where values of 
the K index were as high as 38. The convective area or squall 
line reached maximum intensity at about 0600 GMT while moving south- 
eastward. The area decreased in intensity to shower activity by 
1500 GMT while the cyclone in central Kansas remained stationary. 
There were numerous reports of hail and tornadoes associated with 
this convective system. 

Light rain and rain showers with radar tops below 7.6 km 
(25,000 ft) occurred along and above the warm front associated with 
the first short wave as it moved eastward as a stable wave to a 
position about 250 km off the coast of Maine at the end of the ex- 
periment. 

The cyclone in central Kansas began to intensify and move south- 
eastward into northeast Oklahoma with the eastward movement of the 
second short wave that was located in the Dakotas by 2100 GMT on 24 
April (Fig, 8) . The second major squall line began to develop at 
this time in eastern Kansas within the intensifying low and south- 
westward along the trailing cold front. Values of the K index were 
as high as 42 in the area, indicating extreme instability. The 
squall line stretched from central Illinois to central Texas with 
maximum radar tops reaching 17.1 km (56,000 ft) at 0000 GMT on 25 
April. Severe weather watches were issued as the squall line moved 








(e) 1800 GMT, 24 April 



(f) 2100 GMT, 24 April 


Fig. 7. 


(Continued) 





(a) Surface. 



Fig. 8- 


Synoptic charts for 2100 GMT, 24 April 1975 
(Fucik and Turner, 1975) . 
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southeastward . 

Between 2100 GMT and 1200 GMT (Fig. 9) on 25 April the cyclone 
moved from Kansas eastward into Kentucky with the eastward movement 
of the second short wave perturbation while the associated cold 
front extended southwestward into Texas. The severe squall line 
continued to grow in intensity and length with maximvun radar tops 
exceeding 19.8 km (65,000 ft) and, at 1200 GMT on 25 April, extended 
along an arc almost 1000 km long from eastern Kentucky southwest- 
ward into northern Alabama and westward into Arkansas (Fig, 7) , 

Both hail and tornadoes were reported with the squall line. 

The strong squall line moving southeastward into the south- 
eastern United States was the major weather feature at the end of 
the experiment at 1200 GMT on 25 April, The second trough ran 
north- south through the Northern Plains. The surface anticyclone 
off the Carolines had moved farther eastward out into the Atlantic, 
and cP air was replacing mT air in the Ohio Valley and the south- 
east United States with the movement of the cold front. 
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5. RESULTS 

A "telescoping" approach will be taken in the presentation of 
results. First, the energy budget for the entire experiment area 
will be examined. Next, synoptic-scale energy budgets of smaller 
volijmes enclosing two squall lines will be presented to show the 
relation of these budgets to the overall budget. Finally, spatial 
maps of selected terms in the budget equations will be examined to 
describe features which are lost in the averaging process. In ad- 
dition to the "telescoping" approach, the relation between MDR 
values and energy parameters will be considered. 

These observational results will verify the hypothesis made 
in the objectives — systematic changes do occur in synoptic-scale 
energy budgets of atmospheric volumes enclosing areas of intense 
convection. 

a. Average energetics of the experiment period 

The energetics of the entire AVE IV area, averaged over all 

nine time periods, indicates rather tranquil conditions, but later 

sections will show that this is deceiving. The average kinetic 

energy budget and average total potential energy budget (Table 

4) were obtained by averaging values of budget terms for 168 cen- 

11 2 

tral grid points, an area of 49.0 x 10 m , over all time periods, 
for a total of 1512 points. 

The kinetic energy budget of the surface- to 700-mb layer is 

characterized by conversion from potential to kinetic energy due 

to cross-contour flow, term -C ; near the ground this process is 

K 

enhanced by frictional effects. Weak horizontal and vertical out- 
flow of kinetic energy occur along with negative dissipation, in- 
dicating transfer of energy from grid to subgrid scales of motions. 
The net effect qf these processes is a local reduction of kinetic 
energy in the surface- to 700-mb layer. 

Kinetic energy is being converted to potential energy in the 
700- to 400-mb layer due to air moving across height contours from 
low to high values. Horizontal and vertical outflow occur in this 
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Table 4, Average energy budget for the AVE IV experiment with comparisons to previous research. 

11 2 

The AVE area is 49.0 x 10 m . 


Pressure Layer 
mb 

K 

lO^J 

IL 

at 

w m**^ 

HF 

K 

W m”^ 

VF 

K 

W m"^ 

-=K 
W m”^ 

D 

W 

G 

W m"^ 

"k 

W 

PW 

W m"^ 

HF 

n 

W m~^ 

VF 

n 

W m"^ 

Si 

W m“^ 

200-100 

4.4 

-1.2 

0.1 

-0.5 

-3.0 

1.4 

-30.4 

3.0 

150.4 

649.8 

-123.2 

17.1 

300-200 

5.4 

-0.5 

1.8 

-1.7 

-2.2 

1.8 

4.3 

2.2 

48.9 

1037.9 

-866.9 

-58.7 

400-300 

3.4 

-0.3 

-0.4 

-0.1 

-1.5 

0.8 

25.8 

1.5 

3.3 

427.2 

-415.4 

-103.7 

500-400 

2.2 

-0.3 

-0.3 

0.2 

-0.8 

0.5 

27.2 

0.8 

38.6 

201.2 

-66.1 

-92.2 

600-500 

1.5 

-0.4 

0.1 

0.3 

-0.4 

0.5 

13.1 

0.4 

32.8 

173.0 

-58.0 

-78.0 

700-600 

1.1 

-0.4 

0.3 

0.5 

0.2 

0.2 

-1.8 

-0.2 

8.8 

58.9 

-28.4 

-68.8 

800-700 

0.8 

-0.3 

0.1 

0.4 

0.1 

0.1 

-9.4 

-0.1 

-31.6 

-540.8 

429.9 

-54.5 

900-800 

0.7 

-0.4 

-0.1 

0.4 

1.4 

-1.5 

-24.7 

-1.4 

-43.8 

-787.7 

634.0 

-26.2 

sfc-900 

0.4 

-0.2 

0.0 

0.2 

2.2 

-2.2 

-21.0 

-2.2 

-107.8 

-525.9 

147.8 

-7.4 

Vertical total 

19.8 

-4.0 

1.6 

-0.1 

-4.0 

1.5 

-16.9 

4.0 

99,6 

693.6 

-346.4 

-472.3 

Mature cyclone 

vicinity, surface to 
tropopause, 

(Kung and Smith, 1974) . 

North America during 
short wave 

14.0 

-24.4 

0.5 

27.4 

-37.3 







development, surface 
to lOOmb, (Ward and 
Smith, 1976) . 

Anticyclone vicinity, 
surface to 200mb, 

19.6 

0.5 

-5.2 

0.0 

-4.4 

-0.3 







(Smith, 1973a) , 

- 

-12.3 

12.7 

0.0 

0.4 

0.0 

- 

- 

- 

- 

- 

- 

Mature cyclone 

vicinity, surface 
to lOOmb, 

(West, 1973) . 

■ 


7.3 

-2.5 

17.4 


40.8 

2.5 

227.3 

-1261.6 

460.0 



! 
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layer as in the lower troposphere, but in this layer kinetic energy 
is being transferred from subgrid scales of motion to larger scales. 
A local decrease in kinetic energy is the overall result. 

The kinetic energy budget near the jet stream, the 400- to 
100-mb layer, is characterized by maximum conversion to potential 
energy; term -C equals -3.0 W m in the layer from 200 to 100 mb. 

Ix 

Maximum transfer of kinetic energy from subgrid to grid scales of 
motion occurs in the 300- to 200-mb layer. Vertical inflow and 
horizontal outflow of kinetic energy is occurring with both pro- 
cesses reaching their greatest magnitudes in the 300- to 200-mb 
layer. These combined processes again lead to a local decrease in 
kinetic energy through the layer. 

The kinetic energy budget of the total volume (surface to 100 
rob) can be characterized by conversion from kinetic to potential 
energy, outward horizontal transport of kinetic energy, inward 
vertical transport, the creation of synoptic-scale kinetic energy 
from subgrid-scale motion, and a local decrease of kinetic energy 
with time. 

Table 4 compares the kinetic energy budget for the AVE IV 
experiment with the budgets of previous investigators. The average 
kinetic energy budget for the AVE IV experiment is quite unlike 
budgets obtained for cases of cyclones. For example, the vertical 

total of term -C in the vicinity of an average mature cyclone is 

-2 ^ -2 
27-4 W m (Kung and Smith, 1974) compared to -4.0 W m for this 

study. Results of Smith's (1973a) study of an anticyclone are 

more similar to the present results. Smith reported values of 
-2 -2 

0.4 W m for term -C and 12.7 W m for term HF in the surface- 

K K 

to 200-rab layer. The budget for North America during a period of 
short wave development (Ward and Smith, 1976) is also more similar. 

Positive values of dissipation often are viewed with skep- 
ticism, and the computational uncertainties associated with this 
residual term always should be acknowledged, but several reports 

of positive values occur in the literature. Vincent and Chang (1975) 

-2 

reported a vertical total of 27.4 W m during the developing stage 

12 2 

of an extratropical cyclone for an area about 1.15 x 10 m usxng 
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moving coordinates. Smith (1973b) gave a value of 1.0 W m for the 
surface- to 200-mb layer over North America prior to cyclogenesis; 
he noted that since smaller waves had a tendency to dampen out during 
the period, perhaps their kinetic energy was contributed to the 
larger scales of motion. Kornegay and Vincent (1976) gave a ver- 
tical total of 1.9 W m when a tropical storm interacted with an 
extratropical frontal system. 

While positive values of dissipation are often attributed to 
errors in the input data, Vincent and Chang (1975) and Kornegay and 
Vincent (1976) concluded that realistic errors alone would not pro- 
duce positive values of dissipation for their studies. They sug- 
gested that subgrid- scale energy sources that are not completely 
detectable using synoptic -scale energy data, such as convection 
or eddy momentum transfer, require compensating mechanisms for 
dissipation which if detectable or capable of parametrization would 
reduce or eliminate positive values. It is clear from the synoptic 
discussion that convection and associated turbulence and wave 
phenomena, which cannot be completely resolved in time or space 
with the available data, are widespread and intense during this 
period of study. The topic of dissipation is discussed further in 
following sections. 

Terms of the total potential energy budget during the AVE IV 
experiment are given in Table 4. Since average vertical motion • 
over the area is upward, values of are negative throughout most 
of the atmosphere, indicating conversion from potential to kinetic 
energy. The value of |^jrj/|c^j ^or the surface- to 100-nib layer is 
about 1% which indicates that the limited region exhibits major in- 
teraction with the surrounding atmosphere. The value of 1% is 
similar to values obtained by other investigators (Gommel, 1973; 

Smith and Horn, 1969) . 

The lower troposphere, from the surface to 700 mb, contains 
vertical outflow and horizontal inflow of total potential energy 
while the 700- to 100-mb layer is exporting energy horizontally and 
importing energy vertically. The values of diabatic effects (term G) 
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show a small range throughout the vertical profile. Diabatic cooling 
occurs in the lower troposphere, but generation of potential energy 
due to diabatic heating occurs in the 600- to 200-mb layer. 

The literature contains few studies on total potential energy 
which can be compared with the current study, but West’s (1973) in- 
vestigation of cyclogenesis over the eastern United States is a 
notable exception (Table 4) . His values of terms G, 
were larger than those presently observed, and shapes of vertical 
profiles of the terms were considerably different. 

Based on the previous results, the atmosphere may appear tran- 
quil during the AVE IV experiment, but significant energy processes 
in both time and space have been lost due to averaging. These im- 
portant' events are described in later sections. 

Vertical motion and squall line development 

A first indication of systematic changes in meteorological 
parameters that relate to areas of intense convection can be found 
in fields of vertical motion. Vertical motion fields at 500 mb 
for the nine AVE IV time periods, with superimposed areas of MDR 
values greater than or equal to 4, are given in Fig. 10. More de- 
tailed maps of MDR values are found in Fig. 7 (p, 30) , The first 
major squall line is already well developed and stretches from 
central Illinois to western Nebraska at 0000 GMT on 24 April, but 
the area of maximum upward vertical motion is small and located 
slightly to the rear of the convective area. The area of upward 
motion associated with the surface low and short wave can be seen 
over the northeastern United States. The area of upward motion 
increases greatly in size and intensity by 0600 GMT when the squall 
line is at its maximum intensity. Values as great as -15.0 yb s ^ 
are located over central Missouri, immediately behind the area of 
strongest convection. The area of upward motion at this time is 
the dominant feature in the e3q>erimental area. 

The squall line decreases in areal coverage and intensity 
after 0600 GMT, The most intense convection is occurring over 
southern Missouri and northern Arkansas by 1200 GMT while the area 
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of strongest upward motion^ reduced to -12.0 yb s ^ at 500 mb, is 
located over central Missouri. Downward vertical motion is associ- 
ated with the remnants of the squall line over Missouri and Arkansas 
at 1500 GMT. Positive vertical motion along the Ohio River Valley 
and New England States at 1500 GMT is associated with rainshower 
activity in those areas with MDR values less than 4. 

When the second major squall line develops over eastern Kansas 
between 1800 GMT and 2100 GMT on 24 April, part of the area is ex- 
periencing slight upward motion less than -4.0 yb s areas of 

downward motion also are in the vicinity# but are weak as well. 

An organized squall line, stretching from Illinois to Texas, 
forms by 0000 GMT on 25 April. Increases in upward vertical motion 
to values of -6.5 yb s ^ accompany this development. Peak inten- 
sity of the second squall line occurs at 0600 GMT at which time 
upward vertical motion at 500 mb reaches a value of -16.5 yb s ^ 
over the southern tip of Illinois. Upward vertical motion decreases 
in magnitude by 1200 GMT as the squall line decreases in intensity. 

These results indicate that areas of upward vertical motion 
develop and decay as areas of convection develop and decay, with 
vertical motion somewhat lagging the convection. As suggested by 
Aubert (1957), Danard (1964), and Ninomiya (1971 a and b) , diabatic 
heating associated with convection probably is responsible for in- 
ducing the systematic changes in these vertical motion patterns. 
Since systematic changes in vertical motion are associated with 
areas of convection in the AVE IV experiment, such changes should 
also occur in terms C^, VF^^, and VFj^ of the energy budget equations 
because these are functions of vertical motion. Since horizontal 
divergence is related to vertical motion, systematic changes are 

expected in terms HF and HF„ of the energy equations as well. 

K iL 

This first indication that systematic changes should occur in the 
synoptic -scale energy budgets of volumes enclosing areas of con- 
vection will be examined in much greater detail in following sec- 
tions . 
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c. The energy budget of a developing squall line 

Energy budgets of relatively small atmospheric voliames enclosing 
the squall line that formed at 1800 GMT on 24 April are examined 
in this section. The approach is designed to describe better the 
nature of synoptics-scale energy processes occurring in the vicinity 
of intense convection. These budgets at five time periods are com- 
pared with each other and to budgets for much larger areas. The 
formation, mature, and initial dissipation stages of development 
are studied using this squall line case while an additional squall 
line case is used later to investigate further the mature and 
dissipation stages. 

1) Selection of area boundaries 

The relatively small atmospheric volumes were chosen such that 
they just enclosed the squall line based on MDR values and moved 
as the squall line moved. Vertical dimensions were from the sur- 
face to 100 mb while horizontal boundaries were defined using the 
grid points on which numerical calculations were performed (Fig. 

5, p. 20) . The budget for each area and each 50-mb pressure layer 
was computed as the average of all grid points on or within the 
boundary. It should be emphasized that selection of horizontal 
boundaries was based solely on MDR data; spatial fields of the 
budget terms were not considered in the selection process. As 
spatial fields of the terms will reveal, a different boundary lo- 
cation could affect significantly computed averages in some cases. 
Figures 11 and 12 show two examples of the boundary selection pro- 
cedure while Table 5 gives the number of grid points and associ- 
ated areas for the boundaries used in this squall line case. The 
convection in Ohio and Michigan at 0600 GMT on 24 April (Fig. 12) 
was not included in the volume because it was found to have a 
much shorter life than the remainder of the squall line and ap- 
peared to be a separate system. While the goal was to barely en- 
close areas of convection, at least 20 grid points were used at 
each time. This meant that during initial development, the 
horizontal boundary was centered over the much smaller area of 
convection. Volumes used in this study are approximately the 
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Table 5. Areas for which average energy budgets were computed 
for the squall line forming at 1800 OIT on 24 April. 


Time 


Number of points 


Area 


(m^) 


1800 GMT, 24 April 24 
2100 GMT, 24 April 24 
0000 (34T, 25 April 30 
0600 GMT, 25 April 35 
1200 GMT, 25 April 39 


8.8 X 10^^ 
8.8 X 10^^ 

10.6 X 10^^ 
12.2 X 10^^ 
13.4 X 10^^ 


same size as those used by such authors as Kxmg and Smith (1974) , 
Kornegay and Vincent (1976), and Vincent and Chang (1975) and are 
much smaller than those used by numerous other investigators. They 
retain the desirable aspects of averaging and still enable con- 
sideration of a limited volume. 

2) Average energetics in the vicinity of the squall line 
Vertical profiles of energy budget terms are examined in this 
section. Values for several time periods are plotted on the same 
diagram to show the systematic nature of energy changes occurring 
in the vicinity of the squall line. To this author's knowledge, 
this is the first time that such results have been presented in 
the literature for most terms of the equations of total potential 
and kinetic energy. The tabulated energy budget of the squall 
line at peak intensity is examined in a later section. 

i. Diabatic processes . Vertical profiles of generation of 
total potential energy by diabatic processes (term G) show a cyclic 
nature during the life cycle of the squall line (Fig. 13) . Radi- 
ative effects dominate near the surface with diabatic heating oc- 
curring through 2100 GMT and diabatic cooling afterwards. The 
upper and middle troposphere exhibit a large increase in heating 
through 0600 GMT and a decrease in heating at 1200 GMT with maxi- 
mum values occurring at approximately 450 mb. The time of maximum 
heating, 0600 GMT, corresponds to the time of maximum convective 



G, (lO^W m^)/100 mb 

Fig. 13. Vertical profiles of term G for the squall line of 
25 April. 

activity; the 450-mb value at 0600 GMT is approximately three times 
the value observed at 0000 GMT. Latent heat release associated 
with convection is almost surely the cause for the obseirv’ed pro- 
file in this portion of the atmosphere. Cooling at about 200 mb 
at 0600-1200 GMT is probably due to radiational losses above the 
cirrus shields of cumulonimbus clouds; such cooling has been ob- 
served by Anthes and Johnson (1968) in their study of Hurricane 
Hilda. 

ii. Energy conversion terms . Values of A-conversion (term 
C^) exhibit a remarkable systematic nature (Fig. 14) . Since 
specific volume is always positive, the sign of the teimi is re- 
lated directly to the sign of vertical motion. Slight conversion 
from potential to kinetic energy is observed near the siirface and 
near the jet stream level at 1800 GMT when the sqiiall line is 
beginning to form. This situation changes rapidly as the convec- 
tion intensifies to strong conversion from potential to kinetic 
energy at all levels with a maximum at about 400 mb. 

Conversion decreases somewhat by 1200 GMT as the convection 
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Fig. 14. Vertical profiles of term C for the squall line of 
25 April. ^ 

decreases in intensity- Strong similarity between the vertical 

profiles of terms and G in the middle and upper troposphere is 

evident. The shape of the profile of is consistent with that 

suggested in a numerical study by Danard (1964) . 

The shapes of vertical profiles of term -C (Fig. 15) are 

K 

considerably different from those of term C . and values of -C 

A K 

are much smaller than those of C^. The ratio of 

surface- to 100-mb layer remains approximately 1% at each of the 

five times. Positive values of term -C indicate conversion from 

K 

potential energy to kinetic energy due to cross-contour flow. 

Conversion to kinetic energy decreases with storm development 

and becomes conversion to potential energy at 0600 GMT near the 

800- to 600-irib layer. The opposite occurs above 600 mb where 

conversion to kinetic energy increases with storm development, 

reaches a maximum value at 0600 GMT near 250 mb, and decreases at 

1200 GMT. Results of Mclnnis and Kimg (1972) for term -C do not 

K 

appear as cyclic in nature as the present results. 

Variations in term -C can result from changes in the vector 

K 

wind and/or changes in the height gradient. Convection has been 
related to changes in both parameters in a way that helps esq^lain 
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Fig. 15. Vertical profiles of term -C for the squall line of 
25 April. 

the variation observed in this term.. Increases in values of 
geopotential height in the middle and upper troposphere and decreases 
near the ground together with enhanced convergence in lower levels 
and enhanced divergence in upper levels have been associated with 
thiinderstoirm activity (Aubert, 1957; Danard, 1964; Ninomiya 1971a 
and b) . These combined effects would lead to enhanced flow across 
height contours from high to low values near the surface producing 
enhanced conversion to kinetic energy. In the middle troposphere, 
increased cross- contour flow from low to high values leads to en- 
hanced conversion to potential energy. Enhanced conversion to 
kinetic energy would occur in the upper atmosphere due to in- 
creased flow across height contours from high to low values. 

While the effect of tilting pressure systems has not been considered 
in this e^qjlanation, their presence generally should enhance the 
effects described cibove (Danard, 1964) . Even if no height changes 
occur due to convection, results similar to the cibove could result 
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from wind changes associated with convection together with tilting 
pressure systems. 

Kung (1967, 1974) has shown that conversion to kinetic energy 

is slightly greater at 0000 GMT than at 1200 GMT; the cause is not 

fully understood. Kung's results are opposite the present results 

in the upper atmosphere while at 750 mb Kung found only a small 

diurnal effect. It is doubtful that the systematic nature of term 

-C in this study is caused to any great extent by diurnal effects. 

K 

iii. Boundary fluxes of kinetic energy . Significant sys- 
tematic changes occur in terms HF and VF during the life cycle 

K K 

of the squall line. Horizontal outflow in the lower half of the 

atmosphere and strong horizontal inflow in the upper half of the 

atmosphere occur during the initial formation of the storm (Fig. 16) . 

The horizontal inflow aloft is due largely to horizontal advection 

of kinetic energy since the storms form immediately downwind of a 

jet maximum. The profile reverses as the squall line reaches peak 

intensity at 0600 GMT such that weak inflow in the lower half of 

the atmosphere and strong outflow aloft occur. The jet maximum 

had moved to the northeast of the area of convection by this time. 

At 250 mb the range of HF between 1800 GMT and 0600 GMT is 
2 ^ 

44 W m /lOO mb while in the lower troposphere the range is about 

one order of magnitude smaller. Magnitudes of HF decrease by 1200 

GMT as the squall line decreases in intensity. 

Values of VF are near zero as the convective area forms but 
K. 

change to vertical outflow of kinetic energy below about 400 mb and 
vertical inflow above about 400 mb with further development (Fig. 17) . 
Kinetic energy is being transported to higher levels during the times 
of strongest convective activity since upward vertical motion is 
taking place over the area. 

iv. Dissipation of kinetic energy . Vertical profiles of term 
D, given in Fig. 18, indicate that transfer of kinetic energy from 
subgrid to larger scales of motion generally increases with storm 
development below approximately 600 mb while enhanced transfer from 
grid to. subgrid scales of motion occurs above 600 mb. While some 
cyclic nature is apparent, it certainly is not as pronounced as for 
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Fig. 18- Vertical profiles of term D for the squall line of 
25 April. 

most of the previous terms. Vertical totals of term D presented by 
Mclnnis and Kung (1972) also show considerable fluctuations. Since 
turbulence and other subgrid phenomena are large in thunderstorm 
areas, large values of term D of both signs should be expected. The 
diurnal tendency of term D (Kung, 1967) is small and is not thought 
to be important in explaining these results. Satellite pictures of 
this area reveal the presence of wave phenomena in the cloud pat- 
terns. Their presence could be important in explaining the observed 
values of D, but firm conclusions are not possible. 

Time derivative of kinetic energy . The net effect of 
the various physical processes on kinetic energy within the limited 
volume is rather small because, for the most part, the processes 

cancel each other (Fig. 19) . Values below 500 mb are only about 

—2 ““2 9 K 

±2 W m /lOO mb; above 500 mb they are ±6 W m /I 00 mb. Term 

appears slightly more negative as the convection intensifies, but 
the effect is small. 
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Fig. 19. Vertical profiles of term for the squall line of 
25 April. 

vi. Boundary fluxes of total potential energy . The vertical 

profile of term at 1800 GMT indicates horizontal inflow of 

total potential energy in most levels of the atmosphere (Fig. 20) . 

As the squall line develops the layer of inflow deepens while a 

region of outflow, apparent at 2100 GMT, decreases in depth and its 

3 -2 

maximum value increases in magnitude from about 2.0 x 10 W m /lOO 

3 “2 

mb at 2100 GMT to 7.5 x 10 W m /lOO mb at 0600 GMT. Magnitudes 
of the term decrease in value at 1200 GMT as the convection di- 
minishes in intensity. 

Profiles of term VF^ also reveal a cyclic nature (Fig. 21) . 

The top of the outflow layer near the surface rises from about 900 
mb at 1800 GMT to near 500 mb at 0600 GMT. Correspondingly, the 
inflow layer aloft decreases in depth but peak values increase in 
magnitude during this period. Boundary flux terms indicate that 
potential energy in the lower atmosphere is transported horizontally 
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g. 20. Vertical profiles of term HF„ for the squall line of 
25 April. 



VFj^, (lO^W m“^)/100 mb 


Fig. 21 


Vertical profiles of term VF^ for the squall line of 
25 April. 
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into the convective area, transported vertically due to upward ver- 
tical motion, and then exported in the upper atmosphere. 

3) Energy budget at peak intensity 

i. Discussion . The synoptic- scale energy budget for the 

limited volume enclosing the squall line at peak intensity (0600 GMT) 

reveals very large conversion and transport processes. Table 6 

gives this budget for a volume whose horizontal area is 12.2 x 10^^ 

2 

m (Fig. 11, p. 52) about one-fourth the area of the total AVE IV 

e^qperimental network. Diabatic heating, assumed to be the cause of 

-2 

the other energy changes, reaches a value of 569.0 W m in the 

800- to 200-mb layer. The maximum value in a 100-mb layer is 
-2 

222.4 W m which occurs between 500-400 mb. Krishnamurti and 
Moxim (1971) stated that convective heating should be largest in 
the upper troposphere in situations of deep convection with little 
entrainment. Assuming that all of the heating is due to latent 
heat of condensation, the value of 569.0 W m corresponds to an 
average precipitation rate of 0.076 cm h ^ (0.03 in h over the 
area whidh seems reasonable since active convective cells cover 
only a small fraction of the total area. Fankhauser (1969) noted 
that if total heat production for a heating area is compared to the 
total latent heat production from all thunderstorms in that area, 
the orders of magnitude are comparable so that the net effect of 
storm presence is realized. 

Values of terms C, and -C indicate large conversion from 
A K 

potential to kinetic energy while the small values of j 

indicate a great deal of interaction between the relatively small 

volumes enclosing convection and the surrounding volumes- Although 

-2 

the surface- to 100-rab value of K-conversion is 52.2 W m , a cor- 
responding increase in kinetic energy does not occur in the volume 

because of dissipation and boundary fluxes. In fact, the surface- 

-2 

to 100-mb layer sustains a net loss of kinetic energy of 6.2 W m 
-2 

because 43.9 W m of kinetic energy is transported horizontally 

-2 

out of the region while 14.5 W m is transferred to s\±>grid 
scales of motion. 



Table 6 


Average energy budget for the limited area enclosing a squall line 

11 2 

at 0600 GMT on April 25. The area is 12.2 x 10 m (Fig, 11, p 52) . 


Pressure Layer 
mb 

K 

lO^J 

3K 

3t 

W m"^ 

HP^ 

W 

W 

-S; 

W 

D 

W m“^ 

G 

W 

PW 

W n ^ 

HF„ 
W m ^ 

VFjj 

W m ^ 

W m“^ 

200-100 

4.8 

-2.2 

7.5 

-0.1 

4.3 

0.9 

-88.4 

-4.3 

676.9 

2830.4 

-461.5 

65.4 

300-200 

6,0 

0.2 

28.8 

-11.6 

16.9 

0.4 

2.4 

-16.9 

474.8 

7869.9 

-6208.1 

-318.1 

400-300 

4.4 

-0.1 

10.8 

-8.8 

14.3 

-12.4 

137.3 

-14.3 

-260.6 

4183.9 

-5096.1 

-803.4 

500-400 

3.6 

-1.4 

1.3 

1.8 

11.0 

-9.3 

222.4 

-11.0 

-57.2 

510.0 

-701.5 

-797.7 

600-500 

2.6 

-0.9 

-0.9 

6.4 

4.1 

0.6 

136.3 

-4.1 

-85.9 

-1953.4 

1652.8 

-630.9 

700-600 

1.7 

-1.0 

-1.2 

7.8 

0.5 

5.2 

45.3 

-0.5 

94.7 

-3184.5 

3516.3 

-406.7 

800-700 

1.0 

-0.5 

-0.3 

2.6 

-2.5 

4.4 

25.3 

2.6 

150.6 

-2385.2 

2913.3 

-223.9 

900-800 

0.9 

-0.3 

-1.1 

1.1 

1.3 

-1.5 

-96.8 

-1.3 

-97.9 

-2649.3 

2306.1 

-110.0 

sfc-900 

0.4 

-0.1 

-1.0 

0,7 

2.4 

-2.7 

-145.4 

-2.4 

-279.1 

-2719.4 

1740.1 

-28 ;3 

Vertical total 

25.5 

-6.2 

43.9 

-0.1 

52.2 

-14.5 

238.4 

-52,2 

616.3 

2502.4 

-347.6 

-3253.8 



65 


The vertical total of -C is about twice the value given by 

K 

Kung and Smith (1974) (Table 4, p. 41) for the vicinity of mature 
cyclones and suggests that large areas of intense convection may 
affect the general circulation as much as or more than mature 
cyclones. A striking difference between results of this study 
and those involving matiire cyclones (Table 4) is found in the hori- 
zontal boundary flux term of kinetic energy. Net outflow in the 
surface- to 100-mb layer is seen in this study while large net in- 
flow is observed in most mature cyclones; Petterssen and Smebye 
(1971) have found export associated with a mature cyclone, how- 
ever. While net vertical flux of kinetic energy is near zero for 
the average mature cyclone. Smith (1973b) and West (1973) obtained 
vertical outflow near the ground and inflow in the upper atmosphere, 
as observed in this case. Petterssen and Smebye (1971) , however, 
found vertical export aloft, and import into the middle troposphere. 
The vertical total of term D is larger in the vicinity of mature 
cyclones than near the squall line case at peak intensity. Vertical 
profiles of term D appear similar in mature cyclones and in this 
squall line case with transfer to subgrid scales of motion occiir- 
ring near the ground and near the jet stream, while slight transfer 
of energy to larger scales of motion occurs in the middle tropo- 
sphere (Smith, 1973a and b) . 

The budget of total potential energy is dominated by the boun- 
dary flux terms, HF and VF . As stated previously, sample cal- 

“ 3JT 

culations indicated that term — was several orders of magnitude 

o t 

smaller that terms HF^^ and VF^. Values of HF^ and VF^ given by 
West (1973) for individual layers of a mature cyclone are similar 
in both magnitude and sign to the present results. 

Many of the processes observed near the squall line are re- 
flected in the energy budget for the entire AVE area at 0600 GMT 
(Table 7) . Vertical totals of all terms except G have the same 
sign for both areas, and values at individual layers indicate that 
maxima and minima for individual terms often occur near the same 
pressure level. Since values for the limited volume are often an 
order of magnitude greater than corresponding values for the large 



Table 7 


Average budget for the AVE IV experiment area at 0600 GMT on 25 April. 
The area is 49.0 x lO^ m^. 


Pressure Layer 

K 

3K 

3t 



-s< 

D 

G 


PW 




mb 

5 -2 

10 J m ^ 

W 

W m ^ 

W 

W m ^ 

W m"^ 

W 

W m”^ 

W m“^ 

w 

W m ^ 

W m ^ 

200-100 

4.2 

-0.7 

0.7 

-0.1 

-1.1 

1.1 

-20.2 

1.1 

216.4 

■ 776.4 

-18.7 

22.6 

300-200 

5.6 

0.9 

1.6 

-1.7 

3.5 

-2.7 

1.3 

-3.5 

72.5 

920.6 

-666.7 

-45.2 

400-300 

3.4 

0.1 

0.8 

-1.2 

0.3 

-0.7 

16.0 

-0.3 

88.3 

655.6 

-346.6 

-73.7 

500-400 

2.1 

-0.9 

0.5 

-0.3 

-0.5 

-0.2 

35.2 

0.5 

59.3 

349.4 

-139.9 

-75.4 

600-500 

1.3 

-0.9 

0.3 

1.0 

-0.7 

1,2 

6.5 

0.7 

59.6 

-7.4 

216.2 

-58.6 

700-600 

0.9 

-0.7 

0.5 

1.4 

-0.6 

1.7 

-13.8 

0.6 

65.0 

125.7 

101.8 

-40.5 

800-700 

0.6 

-0.3 

0.3 

0,5 

-0.7 

1.2 

-12.0 

0.7 

7.3 

-288.6 

314.2 

-30.2 

900-800 

0.5 

-0.1 

0.1 

0.3 

0.9 

-0.6 

-81.6 

-0.9 

-6.9 

-329.9 

305.6 

-12.7 

sfc-900 

0.3 

-0.0 

0.1 

0.1 

1.6 

-1.4 

-109.4 

-1.6 

-72.9 

-288.8 

33.0 

-4.0 

Vertical total 

18.9 

-2.5 

4.7 

-0.0 

2.7 

-0.6 

-178.1 

-2.7 

489.3 

1913.1 

-201.1 

-317.7 
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volume, this suggests further the great impact that relatively small 
volumes enclosing convection have on surrounding larger volumes. 

The importance of convection also can be seen by comparing the 
budget of the limited volume enclosing the squall line (Table 6) 
with the budget obtained by averaging over the entire AVE IV area 
and over all nine time periods (Table 4, p. 41). Differences in 
sign of vertical totals of K-conversion and dissipation are especial- 
ly noticeable; vertical profiles of these two terms are quite dif- 
ferent for the two averages as well. Although the large area is 
characterized by K-conversion to potential energy and transfer of 
kinetic energy to grid scales of motion, significant areas of the 
opposite processes are occurring. The remaining terms often ex- 
hibit similar shaped vertical profiles, but the magnitudes of the 
terms for the individual squall line are generally much larger than 
those obtained by averaging over the entire experiment. 

These results indicate that imbedded within the tranquil- 
looking energy budget of the entire AVE IV experiment is a squall 
line whose energy processes are larger than those reported for 
mature cyclones. 

ii. Energy balance of the total potential energy equation . 

When the total potential energy budget is computed for such a 
limited volume and for only one time period, the question of a bal- 
anced budget naturally arises. Good balance is found for this 
study, however. A residual was computed from the potential energy 

0H 

equation, (21); term was neglected because of its small size. 

The residual represents errors due to the input data, the assump- 
tions made, and limitations of the computational procedxires/ and 
also reflects an inability to resolve relatively small-scale 
features in both time and space with the available data. To some 
extent then, it indicates a transfer of total potential energy be- 
tween subgrid scales and larger scales. The value of the residual 
always is smaller than the values of the more important budget 
terms in the nine 100-mb layers and frequently is an order of mag- 
nitude smaller at 0600 GMT on 25 April. The ratio, of the residual 
to the sum of the absolute values of the other individual terms 



68 


in (21) is less than 0.15 in most layers. Balance is generally 
poorest in the 200- to 100-mb layer. A ratio of about 0.12 was 
obtained in 250-rtib layers for an available potential energy study 
of North America during an entire month (Smith and Horn, 1969) . 
Balance for some time periods and volumes of this study is better 
than at 0600 GMT, and worse than at 0600 GMT for others, but it is 
felt that the trends of the results and the nature of the inferred 
processes can be considered with great confidence. 

The energy budget of a decaying squall line 

A second squall line case is examined in this section to 
determine if results of the previous case are an isolated phenomenon. 
This squall line began prior to the start of the AVE IV experiment 
and was already well developed, but not at peak intensity, when the 
experiment began at 0000 GMT on 24 April. The dissipation stage 
of the life cycle of convection can be examined in greater detail 
with this case. The same investigative approach as before was 
used for this squall line; Table 8 gives horizontal areas of the 
limited volumes used to enclose the squall line at four sampling 
times . 

Table 8. Areas for which average .energy budgets were computed 
for the squall line forming at 0000 GMT on 24 April. 


Time Number of points Area (m ) 








11 

0000 

GMT, 

24 April 

22 

00 

to 

X 

10 

0600 

GMT, 

24 April 

28 

10.0 

X 

10^^ 

1200 

GMT, 

24 April 

25 

9.1 

X 

10^^ 

1500 

GMT, 

24 April 

20 

7.6 

X 

H 

r^ 

O 

H 


1) Average energetics in the vicinity of the squall line 
i. Diabatic processes - Heating in the middle and upper 


troposphere is the dominant feature in the vertical profile of 
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-2 

diabatic processes (Fig. 22) . The maximum value of 285 W m /lOO 
mb occurs at 0600 GMT at 450 mb, near the time of peak convective 
activity. Values of heating decrease after 0600 GMT. Radiation 
effects dominate near the ground and near 150 mb in a similar man- 
ner to that observed for the previous squall line case. 

ii. Energy conversion terms . The value of term in the 
surface- to 100-mb layer is near zero at the beginning of the ex- 
periment although some conversion to kinetic energy is evident 
below 600 mb and above 300 mb (Fig. 23) . Conversion to kinetic 
energy increases at all levels at 0600 GMT then steadily decreases 
so that the range of values observed near the time of dissipation 
(1500 GMT) is similar to that observed in the formation stage. 

The maximum value, occurring at about 450 mb, is somewhat larger 
for this case. than for the previous case; the same is true for the 
maximum value of term G. 

Conversion from potential to kinetic energy due to cross- 
contour flow (term -C ) is occurring at all levels at 0000 GMT 
(Fig. 24) . The conversion to kinetic energy increases in magni- 
tude by the time of peak convective activity, then decreases by 
1200 GMT so that conversion to potential energy occurs above 
300 mb. At the last time period (1500 GMT) conversion to kinetic 
energy exists below 400 mb with conversion to potential energy 
occurring above that level. The cyclic nature of changes in 
term -C is similar to that observed for the previous case 
above approximately 600 mb, but the tendency for enhanced conver- 
sion to potential energy in the 800- to 600-mb layer associated 
with peak convective activity is not as evident, however. 

The ratio |”^|r|/|C^j surface- to 100-mb layer is near 

0.35 at 0000 GMT but is less than 0.06 at the remaining times, and 

vertical profiles of -C and C are generally quite dissimilar 

K A 

which indicates major interaction between the limited volume and 
the surrounding volume. 

iii. Boundary fluxes of kinetic energy . Values of term 

HF (Fig. 25) are cyclic in nature and appear very similar to those 

lx 

of the previous squall line. Outflow of kinetic energy in the 



Pressure 



Fig. 23 


Vertical profiles of term C for the squall line 
of 24 April. 




Pressure (mb) Pressure (mb) 



Fig. 24. Vertical profiles of term for the 
squall line of 24 April. 



. Fig 


25 


Vertical profiles of term HF^^ 
the squall line of 24 April. 


for 
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lower half of the atmosphere and inflow in the upper half at 0000 
GMT change to inflow in the lower half and outflow aloft by 0600 
GMT, then return to the original configuration by 1500 GMT. Low- 
level horizontal inflow at peak intensity is more pronounced for 
this case than for the previous squall line. The value of horizon- 
tal transport near 300 mb is influenced by the proximity of the 
convection to the jet stream. Kinetic energy is advected into the 
volume enclosing convective activity during the formation stage 
and advected out of the volume afterwards. 

Values of VF (Fig. 26) also are similar to those observed 
K 

previously. Minimal vertical flux at 0000 GMT changes to strong 
export of kinetic energy in the lower half of the atmosphere and 
strong import of kinetic energy aloft at the time of peak convec- 
tive activity (0600 GMT) . Values then decrease to near zero by 
1500 GMT. Since upward vertical motion occurs in the area, values 

of HF and VF indicate that low-level kinetic energy is imported 
K K 

horizontally, carried aloft, and then exported horizontally near 
the jet stream. 

iv. Dissipation of kinetic energy . Vertical profiles of 
D (Fig. 27) are considerably different for this case than for the 
previous one. Instead of enhanced transfer of kinetic energy from 
siibgrid scales of motion to larger scales near 700 mb, values at 
this level remain nearly constant with no systematic changes oc- 
curring. Above 500 nib, transfer of energy to subgrid scales of 
motion at 0000 GMT gradually changes to transfer of energy to lar- 
ger scales by 1200 GMT, then becomes near zero by 1500 GMT. Phys- 
ical processes that are unresolvable in space and time produce a 
different result for this squall line, but their nature is unknown. 

V. Time derivative of kinetic energy . Values of term 
-2 

range between ±3Wm /lOOmbin the lower half of the atmosphere 

-2 -2 

and are between 3 W m /lOO mb and -12 W m /lOO mb above 500 mb 
* 3k 

(Fig. 28) . In this squall line case, term appears less nega- 
tive when the enclosed convection is intense. The effect is small 
and differs from the result of the previous case indicating that 



Pressure (mb) Pressure (mb) 



-20 -16 -12 -8 -4 0 4 8 12 16 20 

(V7 m ^)/100 inb 

Fig. 26. Vertical profiles of term VP^ for 
the squall line of 24 April. 



Fig. 27. Vertical ^>rofiles of term D for the 
squall line of 24 April. 
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the net balance between the component processes affecting kinetic 
energy is somewhat different. 

vi . Boundary fluxes of total potential energy . Values of 
terms HF^ (Fig. 29) and VF^ (Fig, 30) indicate that total potential 
energy is imported horizontally in the lower half of the atmo- 
sphere, transported aloft, and then exported horizontally when con- 
vection enclosed by the volume is at peak intensity (0600 GMT) . 
Values of the two terms decrease systematically in magnitude after 
that time. The boundary flux processes are similar to those ob- 
served during the previous squall line. 

2) Energy budget at peak intensity 

Table 9 gives a composite energy picture of this squall line 

case at peak intensity, 0600 GMT on 24 April, for an area of 
11 2 

10.0 X 10 m shown in Fig. 12 (p. 53 ) , The magnitudes and types 
of energy processes occurring at this time are similar to those 










Table 9 


Average 
at 0600 


energy budget for the 
GMT on 24 April. The 


limited area enclosing a squall 
area is 10.0 x 10^ (Fig. 12 


line 
p. 53) . 


Pressure Layer 
mb 

K 

10^ J m"^ 

3K 

9t 

W m*"^ 

“■k 

W m”^ 

^■K 

W 

-'=K 
W m“^ 

D 

’"2 

W m 

G 

W m”^ 

W m ^ 

PW 

W 

HFn 
W tn"^ 

W 

W 

200-1CO 

5.2 

-4.3 

10.9 

-3.7 

12.2 

-9,3 

-84.7 

-12,2 

303.9 

2673.0 

-1609.5 

20.4 

300-200 

5.4 

-6.5 

29.5 

-13.7 

11.8 

-2.4 

50.2 

00 

1— 1 

111.3 

7914.1 

-7524.6 

-469.7 

400-300 

3.7 

-2.1 

12.8 

-4.5 

6.2 

0.1 

234.1 

-6.2 

-498.1 

4967.5 

-6710.7 

-1015.3 

500-400 

3.0 

0.2 

3.7 

0.8 

4.8 

-0.1 

285.0 

-4.8 

-501.2 

1187.0 

-2941.4 

-1109.4 

600-500 

2,4 

0.5 

-2.4 

6.,5 

3.3 

1.5 

257.3 

-3.3 

-442.6 

-2683,8 

1134.3 • 

-957.9 

700-600 

1.7 

0.4 

-3.6 

6.8 

4.4 

-0.6 

102.8 

-4.4 

-308.2 

-5536.9 

4457.7 

-683.6 

800-700 

1.3 

-0.6 

-3.2 

3.8 

2,9 

-2.9 

7.6 

-2.9 

-18.3 

-5797.0 

5732.9 

-383.0 

900-800 

1.0 

-1.4 

-2.3 

2.2 

0.5 

-2.1 

-97.8 

-0.5 

-161.2 

-4697.5 

4131.1 

-163.0 

sfc-900 

0.4 

-0.7 

-0.9 

1.0 

1.5 

-2.1 

-105.2 

-1.5 

-192.8 

-2871.6 

2195.4 

-34.1 

Vertical total 

24.1 

-14,4 

44.6 

-0.8 

47.5 

-18.0 

649.4 

-47,5 

-1707.0 

-4845.3 

-1134.8 

-4795.7 
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observed for the previous case at peak intensity (Table 6 , p, 64 ) . 
Heating in the 800- to 200-mb layer totals 937.0 W iti , correspond- 
ing to a rainfall rate of 0.13 cm h ^ (0.05 in h if all of the 
heating is due to latent heat of condensation, and is somewhat 
larger than the previous case. Maximum heating for both cases 
occurs near 450 mb. 

A- conversion to kinetic energy (term C^) is also greater for 
this case than the previous one although the vertical profiles are 
similar in shape. Only about 1.0% of term is realized in cross- 
contour conversion to kinetic energy (term -C ) . The vertical total 

2 ^ 
of 47.5 W m for term -C is much larger than observed in the 

K 

vicinity of mature cyclones (Kung and Smith, 1974) . Vertical 

totals of terms HF and VF are similar for the two cases. Bound- 
K K 

ary flux and dissipation processes are large enough in both cases 

to overcome the large generation of kinetic energy and produce 

3k 

negative values of term in most layers of the atmosphere. 

Boundary fluxes of total potential energy are somewhat larger 
for this case than for the previous one and dominate the budget. 

Some of the processes observed in the relatively small volume 
enclosing the squall line at peak intensity are reflected, to a 
lesser degree, in the budget of the entire AVE area at 0600 GMT 
(Table 10) , but there are some notcible exceptions. The large 
volume is characterized by K-conversion to potential energy and 
transfer of kinetic energy from subgrid to grid scales of motion 
while the opposite processes occur near the area of intense con- 
vection. Diabatic cooling near the surface and above 300 mb is 
large enough to cancel the weak heating occurring in the middle 
troposphere for the larger volume, but heating is indicated in 
the vertical total of the small volume. 

Values of all terms are smaller for the large volume than 
for the more limited volume and like the previous squall line case, 
some terms are an order of magnitude smaller. Vertical profiles 
of the various boundary flux terms and time derivative term are 
generally similar in shape for both volume sizes. 

In summary, most of the systematic changes observed for the 


Table 10 


Average energy budget for the AVE IV es^eriment area at 0600 GMT on 24 April. 

. ^ 2 
The area is 49.0 x 10 m . 


Pressure Layer 
mb 

K 

, -2 
10 J m 

3K 

3t 

W 

W 

w m“^ 

-=K 

W 

D 

W m”^ 

G 

W m”^ 

W 

PW 

W 

W 

W m ^ 

■ 

W m"^ 

200-100 

4.6 

-1.7 

2.7 

0.4' 

-4.5 

5.9 

-66.2 

4.5 

269.4 

943.9 

-1.1 

50.4 

300-200 

5.3 

-1.3 

4.6 

-2.6 

-2.9 

3.6 

-13,9 

2.9 

117.4 

1435.9 

-1025.3 

-45.0 

400-300 

3.3 

-0.1 

0.6 

-0.7 

-1.1 

0.8 

22.4 

1.1 

110.5 

1063.6 

-677.0 

-104.6 

500-400 

2.1 

-0.3 

0.2 

-0.1 

-1.6 

1.5 

20.5 

1.6 

54.3 

570.1 

-380.2 

-120.3 

600-500 

1.5 

-0.4 

-0.1 

0.9 

-2.1 

2.5 

16.0 

2.1 

53.7 

-21.6 

209.5 

-102.0 

700-600 

1.2 

-0.5 

0.0 

0.7 

-1.4 

1.6 

-13.9 

1.4 

51.7 

-259.8 

440.9 

-70.4 

800-700 

1.0 

-0.6 

0.1 

0.3 

-0.9 

0.8 

-39.6 

0.9 

19.2 

-505.7 

572.8 

-42.2 

900-800 

1.1 

-0.3 

0.2 

0.6 

1.2 

-0.7 

-98.6 

-1.2 

-28.5 

-591.5 

491.5 

-16.4 

sfc-900 

0.5 

-0.1 

0.2 

0.3 

2.9 

-2.5 

-100.7 

-2.9 

-82.0 

-203.6 

-84.2 

-2.9 

Vertical total 

20.5 

-5.4 

8.7 

-0.2 

-10.5 

13.5 

-274.0 

10.5 

565.6 

2431.2 

-453.0 

-453.3 
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first squall line case are observed in the second case as well. Al- 

Sk 

though there are consider aible differences in terms and D between 
the two cases, vertical profiles and magnitudes of the remaining 
terms are generally similar. The two sets of results suggest that 
observed energy changes associated with intense convection are not 
isolated phenomena. 

e. Spatial maps of energy budget terms 

The final step of "telescoping in" on energy processes associ- 
ated with convection is a discussion of spatial maps of energy bud- 
get terms. These maps provide a means of directly relating ener- 
getics to observed areas of convection without averaging over par- 
ticular volumes. Since the limited volumes of the previous sec- 
tions were defined solely on the basis of MDR data, average values 
for such volumes do not describe adequately energy processes that 

are not centered over the area of convection. Fields of terms G, 

9k 

C , "C / HF , VF , and D are shown at three time periods on 

A K K K dt 

25 April at 0000 GMT, 0600 GMT, and 1200 GMT, so that energetics 
of the squall line can be examined before, during, and after peak 
intensity. It is not possible to present spatial fields of all 
terms in the energy budget equations at all nine time periods be- 
cause of the large number of maps involved, but those presented 
describe some of the more important energy processes and establish 
time continuity and vertical consistency of the fields. In the 
discussion that follows, layer 1 will be used to refer to the 
surface- to 700-mb layer, layer 2 to the 700- to 400-mb layer, and 
layer 3 to the 400- to 100-mb layer. Positive values of terms are 
shown by solid lines, and negative values by dashed lines; all 

values represent integrated quantities through the given layer 
-2 

e3q>ressed in W m . Convection with MDR values greater than or 
equal to 4 is indicated by scalloping. 

1) Spatial fields at 0000 GMT on April 25 

Diabatic processes . Areas of generation of total poten- 
tial energy due to diabatic heating in layer 2 relate well to the 

precipitation occurring over South Dakota, the Appalachians, and 
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the Midwest States (Fig. 31a) . The convection occurring near 
Oklahoma had foirmed during the previous hour, and no heating is in- 
dicated in that area. Maximum diabatic heating over the squall line 
-2 

is about 600 W m in layer 2 over Missouri. 

ii- Energy conversion terms . A-conversion to kinetic energy 
occurs over many of the precipitation areas in layer 1 (Fig. 31b) 
with maxima located over Missouri and West Virginia. Conversion to 
kinetic energy due to cross-contour flow (term -C ) occurs over al- 
most the entire AVE IV area (Fig. 31c) which is expected since fric- 
tional effects aid in this conversion near the ground. Maximum 

values of term -C are over North Carolina, slightly east of the 
K 

shower activity and southeast of the maximum value of term C^. No 

closed centers of K-conversion are seen at this time over the Mid- 

v/est States associated with the developing squall line. 

Negative values of term C^, indicating conversion to kinetic 

energy, occur in layer 3 along the squall line in the Midwest 

3 -2 

with magnitudes reaching -2.5 x 10 W m over southeastern Okla- 
homa (Fig. 31d) . Areas of maximum K-conversion in layer 3 (Fig. 31e) 
occur near the areas of maximum A-conversion. Conversion to 
kinetic energy is seen along and slightly west of the squall line 
while conversion to potential energy is occurring over the southern 
Mississippi Valley. 

iii. Horizontal boundary flux of kinetic energy . Weak low- 

-2 

level inflow of kinetic energy, less than -10 W m , occurs along 
and just east of the squall line in the Midwest (Fig. 31f) . Areas 
of horizontal inflow and outflow are associated with the older 
precipitation area occurring along the Appalachians. Similar areas 
will be seen along the squall line at later time periods. The con- 
ditions seen in layer 1 generally reverse in layer 3 (Fig. 31g) 

as seen previously in budgets for limited areas. Maximum magni- 

-2 

tudes of term ^ in layer 3 generally are less than 75 W m 
K 

iv. Vertical boundary flux of kinetic energy . Fields of 

VF in layers 1 and 3 (Figs. 31h-i) and vertical motion (Fig. 10, 

K 

p. 45) indicate that kinetic energy is transported aloft near the 

squall line and near the Appalachians. Magnitudes of VF are rather 

K 
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-2 

weak at this time, less than 6 W rti in layer 1 and less than 30.0 
W m ^ in layer 3. 

V. Dissipation of kinetic energy . Most of the area in layer 1 
is experiencing negative dissipation which is expected since layer 1 
includes the surface (Fig. 31j) . Positive values are much smaller 
than the negative values and do not show much relation to the con- 
vection. Centers of positive and negative dissipation are seen in 
layer 3 (Fig. 31k) . Transfer of kinetic energy from siibgrid to grid 
scales occurs in the northeastern portion of the squall line while 
the reverse occurs in the central and southwestern portions. 

vi. Time derivative of kinetic energy . Relatively small 

values of are seen in layer 1; negative values are largest in 
o t 

the Appalachians and are associated with a low-level wind maximum 
(Fig. 311) . Maximum increases in kinetic energy in layer 3 are found 
in northern Illinois just east of an advancing short wave and north 
of the squall line (Fig. 31m ) , but values are relatively small when 
compared to the component energy processes. 

2) Spatial fields at 0600 GMT on April 25 

Although the vertical profiles of energy budget terms shown 
previously indicate considerable changes between 0000 GMT and 0600 
GMT, the changes in spatial fields are even more dramatic. Figure 
32 contains maps for 0600 GMT. 

i. Diabatic processes . The area of diabatic heating in layer 

2 grows in area and intensity as the squall line intensifies (Fig. 

3 -2 

32a). The maximum value of 2.3 x 10 W m is nearly four times 
larger than the maximum value observed 6 h earlier while the area 
of heating corresponds almost perfectly to the area of convection. 

The maximum value corresponds to a rainfall rate of 0.32 cm h ^ 

(0.13 in h ^) which is well below the maximum observed rainfall 
rate of 5.1 cm h ^ (2.0 in h ^) , but as Fankhauser (1969) pointed 
out, such differences are reasonable when area coverage is considered. 
Weak diabatic cooling is occurring over most of the remainder of the 
AVE IV area. 

ii. Energy conversion terms . Maximum A-conversion to kinetic 
energy in layer 1 occurs along and slightly east of the squall line 

















(j) Term D, W m , layer 1 


Fig. 32 . (Continued) 
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(Fig. 32b) . The area has expanded since 0600 GMT, and maximum values 

3 -2 

have nearly doubled to near -1.8 x 10 W m 

A well defined area of K- conversion to kinetic energy (Fig. 32c) 

has formed over western Tennessee, just east of the squall line and 

near the location of maximum A-conversion to kinetic energy. The 

' -2 

maximum value of -C is 44 W m or about 2.5% of the maximum value 

K 

of C^. Relatively weak K-conversion to potential energy occurs over 

much of the area covered by the squall line; centers are located to 

the northeast and southwest of the strong area of conversion to 

kinetic energy. The spatial fields provide valuable information 

in this situation since the maximimi value of term -C is located 

K. 

barely inside the limited area over which averages were computed 
(Fig. 11, p. 52) . 

The area of A-conversion to kinetic energy in layer 3 is much 

larger at 0600 GMT than at 0000 GMT; and maximum values have in- 

3 2 

creased to -3.8 x 10 W m over southwestern Missouri (Fig. 32d) . 

Dual centers are located near the area of most intense convection 

while much weaker conversion to potential energy occurs over most 

of the remainder of the AVE IV area. 

K-conversion to kinetic energy occurs over nearly all of the 

-2 

squall line in layer 3 (Fig. 32e) . The maximimi value of 225 W m 

is found over central Illinois, along the northern edge of the 

squall line. While this value is very large compared to results 

of previous investigators, most previous work has been presented in 

terms of averages only. West (1973) , who has presented spatial 

-2 

fields, reported a maximum value of about 200 W m associated with 
cyclogenesis, however. Synoptic maps for the current study reveal 
wind speeds of 'I'SO m s ^ associated with the jet stream near Il- 
linois with marked cross-contour flow; this suggests that the large 
value is indeed valid. 

iii. Horizontal boundary flux of kinetic energy . Low-level 

horizontal inflow of kinetic energy continues to occur just east of 

the squall line, and maximima values have more than doubled since the 

-2 

last time period and are now -25 W m (Fig. 32f ) . An area of much 
weaker horizontal outflow occurs along the western portion of the 
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squall line. Outside the squall line area, values of HF are rela- 
tively small. 

Very large values of HF„ occur in layer 3 along and northwest of 

the squall line (Fig. 32g) . Maximum values of outflow are near 300 
-2 

W m or about four times the maximum value observed at the previous 

time period. The large magnitude of HF is quite remarkable, and 

^ -4 -1 

values of horizontal velocity divergence, which reach 1.4 x 10 s 

at 200 mb over southern Illinois, are equally remarkable, Ninomiya 

“4 -1 

(1971a) reported divergence values of '^lO s in the vicinity of 
intense convection and found dramatic changes occurring in the value 
as the convection developed. Changes in vertical motion associated 
with the convection have been shown in Fig. 10 (p. 45 ) . Figure 32g 
indicates areas of horizontal inflow of kinetic energy to the 
northeast and southwest of the squall line in layer 3. The presence 
of a jet maximum over northern Illinois helps explain this particular 
configuration of areas. 

iv- Vertical boundary flux of kinetic energy . Vertical export 
of kinetic energy near the ground (Fig. 32h) and vertical import 
of kinetic energy in the upper atmosphere (Fig. 32i) occur at 0600 
GMT near the squall line. Maximum values in both layers have in- 
creased greatly in magnitude while the locations of maximum values 
coincide well with the squall line. Much weaker values occur over 
the remainder of the AVE IV area. 

V. Dissipation of kinetic energy . A pronounced center of 
negative dissipation is present in layer 1 over southcentral Ten- 
nessee, just east of the squall line, while centers of positive 
dissipation occur along the northern and southern portions of the 
squall line (Pig- 32 j). These areas are most likely related to the 
convection (Kornegay and Vincent, 1976) , but an exact explanation is 
difficult, if not impossible, to determine. Ward and Smith (1976) 
stated that downdrafts should add momentum to the grid- scale motions 
near the ground and thus produce positive dissipation, but to conclude 
in this study that positive areas are associated with downdrafts 
and the negative areas with updrafts does not seem warranted. The 
fields of HF and VF in layer 1 (Figs. 32f, h) are not very help- 
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ful in resolving the question of drafts, which is not surprising be* 

cause of the scales of motion that can be resolved. Since fields 

of D are quite similar to those of -C (Fig. 32e) , non-representative 

fields of term -C could be factors in producing the observed fields 
K 

of D. 

Spatial fields of the residual in the vorticity equation have 
been presented by Read (1976) for the AVE IV period. Fields of 
values of his residual are often quite similar to those of term D, 
which also was computed as a residual. The similarity suggests that 
fields of both terms are not caused primarily by computational un- 
certainties, but instead are due mostly to subgrid processes which 
are not resolved in either study since the same data were used in 
each. 

Maximum transfer to subgrid scales occurs over Illinois in 
layer 3 (Fig. 32k) . This area lies to the north of the squall line 
and is flanked by areas of positive dissipation to the east and 
west. Maximum values in these centers are up to three times as 
large as those observed 6 h earlier. The convection undoubtedly 
is related to these areas as well, but again the processes only 
can be speculated. Ward and Smith (1976) suggested that either 
cumulus cells acting as barriers to the flow or nearby jet stream 
activity may cause the observed patterns. 

3k 

vi. Time derivative of kinetic energy . Values of in 
layer 1 (Fig. 321) and layer 3 (Fig. 32m) are small. Increases in 
kinetic energy generally are found east of the squall line while 
decreases are found to the west, but values of the term are often 
an order of magnitude smaller than those of component terms. Although 
small values at 0600 GMT are due partly to use of a 12-h centered 
time difference, trial use of a 6-h uncentered time difference did 
not produce greatly different results. While data are too coarse 
in space to resolve many important features, the problem of temporal 
resolution also should be stressed, even when data at 3- and 6-h 
intervals are used. 

Maximum values of nearly all the terms presented occur near 
the squall line. Clearly, the squall line is the center of major 
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energetic activity in the AVE IV area. 

3) Spatial fields at 1200 GMT on April 25 

The very large values of energy terms observed when the squall 
line was at peak intensity become smaller at 1200 GMT (Fig. 33) as 
the convection decreases in strength. The strongest convection is 
now occurring in southeast Arkansas. 

i. Diabatic processes . Although maximum values of heating 
are approximately half the maximxim magnitude observed at 0600 GMT, 
areas of heating relate well to the location of the squall line 
(Fig. 33a) . Areas of cooling are located to the northeast and west 
of the convective activity. 

Energy conversion terms . Maximum A-conversion to kinetic 
energy continues along and slightly east of the squall line in 
layer 1 (Fig. 33b) while the area of largest values has decreased 
somewhat. An area of conversion to potential energy is located 
to the west of the convective area and has become stronger since 
0600 GMT. 

Maximum K-conversion to kinetic energy in layer 1 is located 
over northern Alabama with conversion to potential energy occurring 
over eastern Tennessee and along the border of Arkansas and Louisi- 
ana (Fig. 33c) . The area over Alabama is smaller and less intense, 
but the areas of conversion to potential energy are defined better 
at this time period. 

The area of maximum A-conversion in layer 3 occurs along the 

Ohio River Valley and over southeast Arkansas (Fig. 33d) . The 

area is not as pronounced at this time period, and maximum values 
3 -2 

are about 25 x 10 W m , about half those observed at 0600 GMT. The 
area of conversion to potential energy continues near Ohio while 
an additional area has developed over Missouri. 

An area of K-conversion to kinetic energy in layer 3 is located 
along the border of Tennessee and Missouri (Fig. 33e) . The area 
has remained nearly stationary and is located to the northwest of 

-2 

the squall line; maximiam values have decreased somewhat to 170 W m 
A short-wave ridge over eastern Texas is related to the area of 
conversion to potential energy located over that area. 



(a) Term G, 10 W m , layer 2 



Fig. 33. 


Spatial fields of energy budget terms at 1200 GMT on 
25 April. Superimposed are areas of MDR values > 4. 
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(m) Term W m ^ layer 3 


Fig. 33. (Continued) 
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iii - Horizontal boundary flux of kinetic energy . Low-level 
horizontal inflow of kinetic energy continues along the eastern por- 
tion of the squall line (Fig. 33f ) . The area of largest values has 
decreased while the greatest observed value has increased. The 
area of weak horizontal outflow to the west of the squall line at 
0600 GMT has increased in intensity and is defined better at 1200 
GMT than at 0600 GMT. 

Values of HF in layer 3 have decreased dramatically (Fig. 33g) . 

^ -2 
Maximum values of outflow are reduced to near 200 W m near the 

-2 

Great Lakes and reduced to near 75 W m over southern Illinois. 
Horizontal outflow occurs generally over the decaying squall line. 

iv. Vertical boundary flux of kinetic energy . Maximum vertical 
outflow of kinetic energy in layer 1 is confined to a relatively 
small area along the Appalachians while an area of vertical inflow 
is located over the central Mississippi River Valley (Fig. 33h) . 

The area of inflow formed during the past 6 h. Opposite conditions 
are found in layer 3 (Fig. 33i) where vertical inflow is located 
near the Appalachians and vertical outflow near Missouri. Vertical 
motion fields (Fig. 10., p. 45) indicate that kinetic energy is 
transported aloft over the Appalachians and downward over Missouri. 

V. Dissipation of kinetic energy . A complex pattern of dis- 
sipation areas is seen in layer 1 (Fig. 33j) . The pattern is some- 
what similar to the previous time period since negative values occur 
ahead of the squall line while positive values are seen along the 
line itself. Positive values have doubled since 0600 GMT so that 
subgrid processes are now a significant source of kinetic energy. 

The patterns and magnitudes of layer 3 near the squall line (Fig. 

33k) are similar to those observed at 0600 GMT; while the squall 
line has moved to the southeast, the centers of term D generally have 
moved east or northeast . 

9k 

vi. Time derivative of kinetic energy . Fields of term 
also are similar to those observed at 0600 GMT in layers 1 and 3, 
but magnitudes have increased (Fig. 331-m) . The areas in layer 1 
have moved little although the squall line has advanced southeast; 
in layer 3, the maximum area over Illinois appears to have moved 
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westward. The shift from a centered to a backward time difference 
is a factor in these motions. 

The decrease in magnitude of most of the energy terms at 1200 
GMT is the most noticeable feature of the previous maps. The forma* 
tion of areas near the squall line where kinetic energy is transported 
downward instead of upward and other areas where conversion to poten- 
tial energy occurs is an interesting phenomenon. Perhaps, the 
atmosphere is undergoing some type of energy readjustment process 
following the violent processes associated with the intense squall 
line. 

The continuity and vertical consistency of the spatial maps 
are further evidence of the reliability of the computations per- 
formed in this study. The maps reveal many interesting processes 
which are lost when averages are obtained, but only partial ejq)lana- 
tions are possible for many processes. 

f • Energy budgets versus MDR values 

1 ) Procedure 

The limited volumes used previously to describe systematic 
energy changes associated with convection contained at least twenty 
grid points where precipitation intensity could vary from MDR 1 to 
9. The relative intensity of convection within each limited area 
was sxibjectively determined from MDR data within those areas. The 
purpose of this section is to relate energetics objectively to the 
intensity of convection by averaging energy parameters over grid 
points that have equal MDR values. This is not a traditional spatial 
average since grid points comprising the average for a particular 
intensity are not necessarily adjacent. MDR values were assigned 
to each numerical grid point (Fig. 5, p. 20) at the nine time pe- 
tiods by taking the maximum value within 1/2 grid distance (^80 km) . 
Averages were then computed over five MDR categories, combining all 
time periods, where MDR 0-9 included all 1512 grid points as de- 
scribed previously, MDR 0-1 included 1114 grid points associated with 
non convection, MDR 2-9 consisted of 398 grid points associated with 
all categories of convection, MDR 4-9 consisted of 193 points 
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associated with thunderstorm and severe thunderstorm activity, and 
MDR 8-9 consisted of 65 grid points associated with severe thunder- 
storm activity. The overlapping classification scheme was used so 
that averages were based on a large number of individual grid points. 
2) Terms of the budget equations 

Results of this section will demonstrate that magnitudes of 
energy processes are related to intensity of convection deter- 
mined from individual grid points. 

i. Diabatic processes . Diabatic processes associated with 

the average of all grid points and grid points of nonconvective 

areas are much smaller at all levels than those observed for grid 

points of convective areas (Fig. 34) . Heating in the middle and 

upper troposphere is pronounced and increases with the intensity of 

-2 

convection reaching a maximum of 140 W m in the 500- to 400-mb 
layer. Since more intense convection occurs in the late afternoon 
and evening, radiative cooling is very large near the surface for 
these categories. Radiative cooling from the tops of cirroform 
clouds occurs in the 200- to lOO-rtib layer . 

ii. Energy conversion terms . Net A-conversion of energy is 

near zero in the layer from the surface to 100 'mb for grid points 

in nonconvective regions, but conversion to kinetic energy steadily 

increases as averages focus in on more intense convection (Fig. 35) . 

The profile of term -C indicates conversion to kinetic energy near 

K 

the ground and a near zero value in the 800- to 700-mb layer for 

all categories (Fig. 36) . Conversion to potential energy becomes 

greater with altitude above 700 mb for the total area and the points 

associated with nonconvective regions, but conversion to kinetic 

energy increases above 700 mb for the grid points associated with 

convection. Maximum values of -C occur near the jet stream in the 

K 

case of MDR values 2 to 9, but reach a maximum near 450 mb for MDR 
values 4 to 9. A doiable maximum, occurring at 450 mb and 150 mb, 
is evident for the category of most intense convection and was not 
observed in previous results. Conversion to potential energy near 
700 mb that was observed in the limited volume enclosing the squall 






Pressure (mb) 


111 



-Cj^r (W m"^)/100 nib 


Fig. 36- Vertical profiles of term -C^ for various ranges of MDR 
values. 

line of 25 April at peak intensity (Fig. 15, p. 57) is not as 
pronounced in this averaging procedure. 

iii. Boundary fluxes of kinetic energy . Values of term HF^ 
show relatively little variation among the five categories below 
700 mb although the category of most intense convection is associ- 
ated with slightly larger inflow (Fig- 37) . Values in the 300- 
to 200-mb layer indicate a systematic increase in horizontal out- 
flow with an increase in the intensity of convection. The level 
separating inflow from outflow is lower for the strongest convec- 
tion than for the average of points having MDR values 2 to 9- This 
feature is not evident in the averages for particular volumes 
that were presented earlier. 

Only small vertical flux is evident in the average of the en- 
tire AVE IV area and for grid points associated with nonconvec- 
tive regions (Fig. 38) , Low-level vertical outflow and upper- 
level inflow increase in magnitude for categories of more intense 
convection. 
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iv. Dissipation of kinetic energy > Transfer of kinetic energy 
to subgrid scales of motion occurs below 800 mb in all five categories 
(Fig. 39) . Above 800 mb, transfer to larger scales of motion occurs 
over the entire area and for grid points associated with noncon- 
vective regions. The profiles for the categories of convection 
show maximum transfer of energy to the larger scales at about 650 
mb and transfer to subgrid scales above 550 mb. Differences between 
categories of convection are pronounced near 650 mb and above 300 mb. 
Since vertical profiles shown previously for limited volumes en- 
closing the squall lines did not yield consistent results between 
the two cases for term D, the present results, which are similar 
to those observed for the squall line of 25 April (Fig. 18, p. 60), 
are encouraging- Although physical processes producing the observed 
profiles of term D are not completely understood, subgrid processes 
of all types are large in areas of convection- Several hypotheses 
to explain the values have been given previously. 

V. Time derivatives of kinetic energy . All five categories 
show losses of kinetic energy (Fig. 40) , but the most intense con- 
vection is characterized by considerably larger losses above 500 mb 

3 K 

than the others. Values of are nearly zero at 750 mb for MDR 

o t 

values from 8 to 9 while they are more negative for the other cate- 
gories of convection. 

Another approach to the time derivative question is to com- 
pute the total derivative which indicates changes following a par- 
cel instead of the local derivative which indicates changes at a 

fixed point- The data in this section are used to compute an average 

DK 

"synoptic- scale total derivative" of kinetic energy, for each 

MDR category- The total derivative can be expressed symbolically as 


DK 

Dt 


-c^ + D, 


(26) 


which indicates that changes in kinetic energy of a "synoptic- scale 
parcel" are due to cross-contour conversion of energy and dissipa- 
tion of energy- Smith (1973b) called the Siam -C + D the net in- 

K 

ternal source or sink of kinetic energy. Figure 41 reveals rela- 

DK 

tively small values of •— for the entire AVE area and for grid 
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Table 11. 


Average energy budget for all grid points with MDR values 4 to 9. 
The number of grid points is 193. 


Pressure Layer 
mb 

K 

-2 

10 J m 

3k 

3t 

W 

HF^ 
W m ^ 

VF^ 

W 

W m ^ 

D 

W 

G 

W m"“ 

W m”^ 

PW 

W m”^ 

HF„ 

W 

VF„ 

W 

W m'^ 

200-100 

4.8 

-2.8 

1.9 

-1.0 

1.1 

-3.0 

CO 

1— t 
1 

-1.1 

279.9 

1369.8 

-390.3 

17.4 

300-200 

5.8 

-2.4 

8.9 

-5.6 

3.1 

-2.2 

10.0 

-3.1 

162.7 

3878.2 

-3308.8 

-189.6 

400-300 

4.2 

-0.2 

3.4 

-3.5 

3.0 

-3.3 

91.4 

-3.0 

-151.9 

2212.3 

-2743.4 

433.4 

500-400 

3.2 

-0.4 

1.5 

0.2 

3.1 

-1.8 

121.0 

-3.1 

-92.5 

593.6 

-917.6 

-451.7 

600-500 

2.6 

-0.6 

0.1 

2.9 

2.7 

-0.2 

80.0 

-2.7 

-102.0 

-805.8 

440.5 

-383,8 

700-600 

1.8 

-0.8 

0.1 

3.5 

2.0 

0.8 

23.8 

-2.0 

-12.4 

-1523.5 

1480.0 

-276.5 

800-700 

1.2 

-0.6 

0.0 

1.9 

0.6 

0.6 

0.2 

-0.6 

-7.7 

-1649.5 

1622.3 

-179.9 

900-800 

0.9 

-1.1 

-0.6 

1.0 

1.6 

-2.3 

-76.9 

-1.6 

-178.1 

-2546.8 

1922.0 

-97.4 

sfc-900 

0.4 

-0.6 

-0.5 

0.5 

2.2 

-2.8 

-104.2 

-2.2 

-308.1 

-2442.5 

1361.5 

-24.3 

Vertical total 

24.9 

-9.5 

14.8 

-0.2 

19.3 

-14.2 

103.8 

-19.3 

-410.2 

-914.3 

-525.8 

-2019.3 



Table 12. Average energy budget for all grid points with MDR values 8 to 9. 
The member of grid points is 65. 


Pressure Layer 
mb 

K 

lO^J m'^ 

3K 

3t 

W m"^ 

HF^ 
W m”^ 

W 

W m"^ 

D 

W m"^ 

G 

W 

W m“^ 

PW 

W 

W m’^ 

W m"^ 

W 

200-100 

5.4 

-5.1 

5.7 

-1.9 

7.8 

-9.0 

“70.1 

-7.8 

478.9 

2239.8 

-563.3 

33.1 

300-200 

6.6 

-7.5 

14.0 

-7.2 

4.2 

-4.8 

11.2 

-4.2 

231.5 

5073.9 

-4263.9 

-259.5 

400-300 

4.5 

-2.9 

8.6 

-3.6 

5.0 

-2.9 

111.5 

-5.0 

-89.7 

3140.1 

-3454.0 

-553.7 

500-400 

3.2 

-1.7 

5.8 

0.3 

6.4 

-1.9 

137.1 

-6.4 

-45.2 

1549.2 

-1707.4 

-592.7 

600-500 

2.4 

-0.8 

2.5 

2.8 

4.4 

0.1 

88.3 

-4.4 

-123.8 

-435.4 

2.0 

-536.5 

700-600 

1.6 

-0.3 

0.7 

5.0 

2.3 

■ 3.2 

35.8 

-2.3 

-4.6 

-1623.6 

1607.4 

-410.8 

800-700 

0.9 

-0.1 

-0.4 

2.3 

0.2 

1.7 

3.6 

-0.2 

-23.9 

-2247.8 

2163.8 

-284.9 

900-800 

0.8 

-0.8 

-1.8 

1.4 

1.0 

-2.2 

-123.1 

-1.0 

-315.1 

-4160.8 

3055.2 

-162.0 

sfc-900 

0.3 

-0.5 

-1.0 

0.7 

1.5 

-2.3 

-160.7 

-1.5 

-374.4 

-3784.2 

2470.4 

-37.3 

Vertical total 

25.7 

-19.6 

34.2 

-0,2 

32.7 

-18.3 

33.5 

-32.7 

-266.3 

-248.8 

-689.9 

-2804.3 
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to potential energy occurring in the 800- to 300-mb layer. While 

this phenomenon was observed to some extent for the squall line of 

25 April (Table 6, p. 64) , it is not apparent in the MDR averages 

of either category. Maximum conversion to kinetic energy is observed 

near the ground and near the top of the atmosphere in both studies. 

Low- and middle- level horizontal inflow of kinetic energy is 

seen in the results of M&K while a much more shallow layer is seen 

for the current study. In fact, values of BF for MDR 4 to 9 are 

K 

near zero below 500 mb. Maximum horizontal outflow near 300 mb is 

common to both studies, but the value given by M&K is about 1/2 

that observed for MDR values 4 to 9 and about 1/4 that observed for 

-2 

MDR values 8 to 9. M&K reported a vertical total of -8.88 W m 

(inflow) while the value for the most intense convection of this 
-2 

study is 34.2 W m (outflow). 

Profiles of term VF are quite similar between the two studies 

K 

below 400 mb, but are considerably different above that level. The 
differences are probably due to the method in which vertical motion 
was computed and due to dissimilar synoptic conditions. 

Values of term D are comparable between the two studies below 
about 500 mb; this is especially true of the average obtained for 
MDR values 8 to 9, Larger negative dissipation at the surface was re- 
ported by M&K, however. M&K showed transfer of energy from si±>- 
grid to grid scales of motion occurring in the 800- to 300-mb layer 
while this feature extends only to 500 mb for the category of most 
intense convection of this study. Both studies indicate maximum 
transfer of kinetic energy to subgrid scales of motion in the 200- 
to 100-mb layer, but the value of M&K was about three times larger 
than that observed for MDR values 8 to 9. A major conclusion by M&K 
was that areas of thunderstorms can dissipate more energy than is 

created and therefore must import the remainder, but current re- 

9 K 

suits do not support such a conclusion. Values for term were 
not given by M&K; they were stated to be small. 

General agreement exists between the present results and those 
of Mclnnis and Kung although there are interesting differences. No 
complete budgets of the total potential energy or available 
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potential energy of volumes enclosing convective areas occur in 
the literatxire to this author's knowledge. 

Energy budget summary 

Average kinetic and total potential energy budgets (Figs. 44-45, 
respectively) of limited volumes enclosing the squall lines at three 
stages in their life cycles are given in schematic form to summarize 
the complex energy processes described in this report. The limited 
volumes used here are identical to those described previously; their 
dimensions are given in Table 5 (p. 54) and Table 8 (p. 68). Values 
for the development stage are an average of those observed in the 
vicinity of the first squall line case at 1800 GMT on 24 April and 
the second case at 0000 GMT on 24 April while values for the matiire 
stage are an average of those observed at 0600 GMT on 24 April and 
25 April- The average of values for the limited volumes enclosing 
the convection at 1200 GMT on 25 April and 1500 GMT on 24 April is 
used for the dissipation stage. Since the budgets for each indivi- 
dual squall line have already been discussed in great detail, only 
general comments are given in this section. 

During the development stage of the squall line, the limited 
atmospheric volume is characterized by horizontal inflow and ver- 
tical outflow of kinetic energy near the level of the jet stream, 
and horizontal outflow and vertical inflow below that level. K- 

conversion to kinetic energy, transfer of kinetic energy to sub- 

9 K 

grid scales of motion, and negative values of term occur through- 
out the vertical extent of the limited volume. Diabatic cooling 
is indicated below 400 nib while heating is indicated in the 400- to 
100-mb layer. A-conversion to potential energy occurs in the 700- 
to 400-mb layer, but conversion to kinetic energy occurs above 
and below that layer. Boundary fluxes of total potential energy 
are relatively small at this stage of development. 

Magnitudes of most energy processes are much larger during 
the mature stage of the squall line than the development stage. Dia- 
batic heating is pronounced in the middle troposphere and is almost 
certainly the cause of the other changes in energy processes. 
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A- conversion and K- conversion to kinetic energy are found in all layers 

of the atmosphere with the maximum value of occurring in the 700- 

to 400-mb layer and the maximum value of -C occurring in the 400- 

K 

to 100-mb layer. The processes of horizontal inflow and vertical 
outflow of kinetic and total potential energy in the lower levels 
of the atmosphere, together with horizontal outflow and vertical 
inflow aloft carry large quantities of kinetic and total potential 
energy up and out of the convective areas. Transfer of kinetic 
energy to subgrid scales of motion occurs in all layers but is some- 
what smaller during this stage than the development stage. Local 
changes of kinetic energy are small and negative in all layers of 
the atmosphere. 

Most energy processes decrease in magnitude as the squall line 
dissipates. Diabatic cooling is observed near the ground and near- 
the level of the jet stream, but diabatic heating is observed in 
the 700- to 400-mb layer. A-conversion and K-conversion to kinetic 
energy continue at this stage, but magnitudes are smaller than 
during the mature stage in most layers of the atmosphere. -Total 
potential energy still is concentrated horizontally in the lower 
atmosphere, transported aloft, and e:?q)orted near the jet stream, 
but magnitudes of these processes also are smaller. Horizontal ex- 
port of kinetic energy occurs in all layers of the atmosphere during 
the dissipation stage while vertical export of kinetic energy occurs 
below 400 mb and vertical import occurs above that level. Kinetic 
energy continues to be transferred to subgrid scales of motion, but 

magnitudes are slightly smaller at this stage than the mature stage. 

9K 

Values of — are positive in the middle troposphere but are negative 
<3t 

above 400 mb and below 700 mb. 

While random errors in input data and computational uncertain- 
ties are known to exist, it should again be emphasized that the 
trends of the results and the nature of the inferred processes can 
be considered with great confidence. This conclusion is suggested 
by the systematic nature of the results, comparisons with other in- 
. vestigators, and the good balance of the potential energy budget. 
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Random error estimates of this study and those of Vincent and Chang (1975) 
and Kornegay and Vincent (1976) also are useful in this regard. 
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6 . -CONCLUSIONS 


Synoptic-scale budgets of kinetic and total potential energy 
for volumes enclosing intense convection have been studied in this 
report. The availability of 3- and 6-h synoptic-scale data from 
the AVE IV experiment enabled temporal resolution of features that 
would ordinarily be lost. 

Several conclusions are drawn from the present research: 

1) Systematic temporal changes occur in most terms of the 
energy budget equations for relatively small volumes that enclose 
intense convection. These changes are related to the life cycle of 
the enclosed convection . 

2) The kinetic energy budget of a volume enclosing a squall 
line at peak intensity is characterized by large values of K- 
conversion to kinetic energy. The processes of horizontal inflow 
and vertical outflow in the lower levels, together with horizontal 
outflow and vertical inflow aloft transport kinetic energy up and 
out of the convective areas. Transfer of kinetic energy to sub- 
grid scales of motion occurs in the volume, while local changes in 
kinetic energy are relatively small - 

3) The total potential energy budget of a relatively small 
volume enclosing a squall line at peak intensity is characterized 
by large values of latent heat release in the middle troposphere . 
Values of A-conversion are large, but values of conversion due to 
cross-contour flow are much smaller indicating large interaction 
of the relatively small volumes with the surrounding atmosphere. 
Total potential energy is concentrated horizontally in the lower 
levels of the atmosphere, transported aloft, and then exported 
near the jet stream. Only small local changes of total potential 
energy occur , however . 

4) These results and those of other investigators suggest 
strongly that the systematic energy changes are due to latent heat 
release associated with convection. It is impossible, however, to 
determine what energy changes would have occurred under similar 
synoptic conditions but without the convection. 
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5) Some of the synoptic- scale energy processes associated with 
the squall line at peak intensity are as large or larger than those 
observed in the vicinity of mature cyclones • Relatively small areas 
of intense convection/ therefore , can have a large effect on larger 
surrounding atmospheric volumes, 

6) ^>Hiile errors in input data and in computational procedures 

lead to errors in the energy terms / these errors are not large enough 

to change conclusions based on those terms. Average energy budgets 

11 2 

for areas as small as 8.0 x 10 m are reasonable/ and spatial 
fields of dnergy budget terms shov/ good time continuity and vertical 
consistency . 

7) Values of energy processes are related to the intensity of 
convection at individual points. The energy budgets of points as- 
sociated with nonconvection differ greatly from those where con- 
vection is observed, 

8) Results of this investigation are in general agreement 
v;ith those of other studies , but since marked variability likely 
exists betv/een individual areas of intense convection, more research 
of this nature should be conducted. 
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7. RECOJ4MENDATIONS FOR FURTHER RESEARCH 


Further studies of the type presented can contribute to a 
better understanding of the energetics of convective systems and 
their interaction with other scales of motion. Since the time and 
length scales of convective systems are relatively small, special 
observational networks should be established to provide data for 
these experiments. Data from projects SESAME and HIPLEX and future 
AVE experiments should be especially useful. 
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APPENDIX 


Rav;insonde Stations Participating in AVE IV Experiment 


Station 

Number 

Location 

208 

(CHS) 

Charleston, South Carolina 

211 

(TP A) 

Tampa, Florida 

213 

(AYS) 

Way cross , Georgia 

220 

(VPS) 

Apalachicola, Florida 

226 

(CEM) 

Centerville, Alabama 

232 

(BVE) 

Boothville, Louisiana 

235 

(JAN) 

Jackson, Mississippi 

240 

(LCH) 

Lake Charles , Louisiana 

248 

(SHV) 

Shreveport, Louisiana 

255 

(VCT) 

V ic tor ia , Texas 

260 

(SEP) 

Stephenville , Texas 

261 

(DRT) 

Del Rio, Texas 

265 

(MAP) 

Midland, Texas 

304 

(HAT) 

Hatteras, North Carolina 

311 

(AHN) 

Athens , . Georgia 

317 

(GSO) 

Greensboro, North Carolina 

327 

(BNA) 

Nashville, Tennessee 

340 

(LIT) 

Little Rock, Arkansas 

349 

(Ul^l) 

Monette, Missouri 

363 

(AMA) 

Amarillo, Texas 

402 

(VVAL) 

Wallops Island, Virginia 

405 

(lAD) 

Sterling, Virginia (Dulles Airport) 

425 

(HTS) 

Huntington, West Virginia 

429 

(DAY) 

Dayton, Ohio 

433 

(SLO) 

Salem, Illinois 

451 

(DDC) 

Dodge City, Kansas 

456 

(TOP) 

Topeka, Kansas 

486 

(JFK) 

Fort Totten, Nev; York (Kennedy Airport) 

518 

(ALE) 

Albany, Nev/ York 

520 

(PIT) 

Pittsburgh, Pennsylvania 

528 

(BUF) 

Buffalo, New Y'ork 

532 

(PIA) 

Peoria, Illinois 

553 

(OMA) 

Omah a, Nebraska 

562 

(LBP) 

North Platte, Nebraska 

606 

(p:-m) 

Portland, Maine 

637 

(PNT) 

Flint, Michigan 

645 

(GRB) 

Green Bay, Wisconsin 

654 

CHUR) 

Huron, South Dakota 

655 

(STC) 

St. Cloud, Minnesota 

662 

(RAP) 

Rapid City, South Dakota 

11001 

(MFS) 

Marshall Space Flight Center, Alabama 

22002 

(FSI) 

Fort Sill , Oklahoma 
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